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ABSTRACT 


Magnetotelluric measurements were made in the 
period range 1 - 1,000 seconds, at sixteen different 
sites in the plains of central Alberta and in the 
Rocky Mountains. Tensor impedances, tensor 
apparent resistivities in principal directions, 
Magnetic and telluric field polarization parameters, 
and coherencies between the vertical component of 
the magnetic field and horizontal components of 
the magnetotelluric fields were computed for each 
Site using digital spectral analyses techniques. 

The results indicate that the magnetotelluric 
fields in the Plains of central Alberta are 
controlled by the strike of the contacts of the 
undeformed sedimentary rocks of the Plains an 

the one hand with the crystalline rocks of the 
underlying Precambrian basement and the deformated 
sedimentary rocks of the Rocky Mountains on the 
other. One dimensional models derived from the 
apparent resistivities for E-polarization are 
consistent with the known sedimentary section. 
Magnetotelluric results in the Rocky Mountains and 
in the foothills are not invariant and are assumed 


to be affected by local three-dimensional structures. 
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CHAPTER 1 
INTRODUCTION 

1.1 Purpose of the Magnetotelluric Sounding 

Since the inception.of.the-magnetotelluric method 
in 1950, the Province of Alberta has been-an important 
experimental site for this method. Several investigators 
(Garland.-and Webster,.1960;-Srivastava,.et al., 1963; 
Vozoff et. al.,-1963; Rankin and Reddy, 1969; Peeples, 
1969) conducted. magnetotelluric experiments in central 
and soughern Alberta.. A considerable ambiguity exists 
in- the. published results since, in addition to an 
expected. variation from location to location, there 
exists-a considerable variation in the published results 
from data obtained at the same site from one investigation 
to another.. In the last decade considerable advances 
have been made-in. understanding the problems involved 
when.:complex-magnetotelluric fields impinge on complex 
geologic.structures. Together with the improvement in 
statistical. methods, this understanding should lead to 
a better and more consistent interpretation of the sub- 
surface of the earth. A regional magnetotelluric survey 
has been conducted in central Alberta whose prime 
objective is to apply these modern developments to 
investigate the electrical conductivity structure of 


-the Earth's crust in this region. 
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1.2 Historical Review of the Magnetotelluric Method 

Tikhonov (1950), and almost simultaneously, Kato 
and Kikuchi (1950) and Rikitake, (1950, 1951), observing 
the correlation between the electric and magnetic field 
variations at the surface of the earth, considered the 
possibility of using these variations in order to find 
the conductivity at great depth... Cagniard.,(1953) 
showed that the interpretation could be carried out in 
a practical manner regardless of the depth by applying 
the law of similitude. This method consists of matching 
the experimental results, plotted to a suitable scale, 
with a set of master curves. 

The master curves of Cagniard are constructed by 
solving the boundary value problem for a horizontally 
stratified earth under the assumption of plane electro- 
magnetic waves incident on the surface of the earth. It 
is easily shown that» in practice, the results.are 
independent of the angle of incidence. 

Wait (1954) has developed a first order correction 
to the simple Cagniard theory for the effect of a finite 
source, and has concluded that these corrections are 
required for sources which vary over horizontal distances 
ooh thiexonde mio fhvomishess) than caniskim depth.) Price (1962) 
has developed a more rigorous theory for finite sources 


from which corrections can be computed in certain cases. 
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However, Madden and Nelson (1964), considering some 
realistic earth models, have shown that for cases of 
practical interest, the plane wave assumption is valid. 
Srivastava (1965) has shown that the effect of the 
source dimension can be neglected for periods less than 
1000 seconds when determining moderate conductivities 
to shallow depths (10 - 20 km.) in the earth. 

Following Cagniard's paper, magnetotelluric 
measurements have been made by many investigators 
(Cantwell, 1960; Wiese, 1962; Srivastava et al., 1963; 
Vozoff et al., 1963; Tikhonov and Berdichevskii, 1966; 
Hopkpnsieand Smieth polo 66 - (Swaitt, 967+ Adams ¥ P2969: 
VOZOLE let cal Ay oRoGo: Koveun and ‘Chicherina, rl1969; and 
others). In the last few years attempts have been made 
by Wu (1968), Chetaev (1966), Nabetani and Rankin (1969) 
and Patrick and Bostick (1969), to solve the inverse 
problem for direct interpretation of magnetotelluric 
data for one dimensional models. 

The relationship between magnetic and telluric 
fYelds iat cthe isurface: for variGus’ ‘types “of“anisotropy 
Lrmtconductipnity ln. -simiglie or imultijlayérec earth “has 
been studied by Chetaev (1960, Cantwell (1960), 
Pokyszeansk 12. 8196191 Bos oi ch! and? Smith? “(1 962) "Mann C2965) 
O'Brien and Monaeneoth (2967) VePraus*-and- "Petr (M69), 


Rankin and Reddy (1970) and others. These authors 
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expressing surface impedance/admittance as a tensor 
quantity, have provided methods to find the principal 
directions of conductivity anisotropy and the conductivity 
ratio between the major and minor axes of anisotropy in 
the case of anisotropic half-space, or in the case of 
multilayered half-space when only one of the layers is 
anisotropic. Madden and Nelsén (1964), Swift (1967), 
Morrison et al., (1968), Sims and Bostick (1969 have 
presented methods to compute tensor impedance elements 

by uSing statistical and spectral techniques. 

The response of the magnetotelluric field to 
vertical contacts and cylindrical inhomogeneities has 
also been examined in the last decade. Neves (1957) 
calculated apparent resistivities over dipping interfaces 
using finite difference techniques. Kovtun (1961) 
has discussed, ina general way, the magnetotelluric 
fields for two-dimensional inhomogeneous structures. 
d"*Ercevilte~and Ktnetz (T1962) "and Rankin (1962) gave 
the analytical solutions for vertical fault and vertacal 
dyke respectively. In both the cases the fault and 
dyke is of finite extent underlain by a homogeneous 
and isouvropic, half-space, andloniy, the case of H- 
polarization (magnetic: field vector parallel to the 
strike) is considered. Weaver (1963) solved the 


problem of infinite depth vertical contact for E- and 
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H-polarizations, numerically. Blake and Swift (1967) 
extended the work of Rankin (1962) and Weaver (1963) for 
adyke Or va long promontory. These latter results are 
incorrect for E-polarization case because of the 
improper boundary conditions used at the air-earth 
boundary in their calculations. Madden (see Swift, 1967) 
obtained solutions for two-dimensional problems via a 
transmission line analogy and Patrick and Bostick (1969) 
developed a finite difference scheme for two-dimensional 
magnetotelluric models. 

Analog model studies for two and three dimensional 
bodies have been made by Rankin et al., (1965), Dosso 
(1966), Takaes (1969) and others. These authors 
Suggested that the analog model studies can successfully 
be used to interpret magnetotelluric measurements over 


two and three dimensional structures. 


1.3 Magnetotelluric Source Fields 


The geomagnetic pulsations in the frequency range 
1 = .001 Hz are the source fields for the magnetotelluric 
method presented in this thesis. The pulsations in 
this range are classified into two main classes; those 
of regular and continuous appearance called Pe, and 
tnoseeween an “irregu Var Storm called -Pir*+ The Pe group 


is further subdivided into five subgroups and the re 
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group into two subgroups (Jacobs et al., 1964). These 
groups and subgroups are presented in Table 1. The 


geomagnetic pulsation spectrum is shown in Fig. 1.1. 


Table 1 


Classification of Geomagnetic Micropulsations 


Type Period 
(Seconds) 


Continuous pulsations 
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It is generally accepted that continuous pulsations, 
except Pc 1, are generated as resonance oscillations 
in the magnetosphere of the earth under the influence 
Of Solar corpuscular streams. “The appearance or dit terent 
types of Pc depend.on the characteristic of solar 
corpuscular streams which generate them, 
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Figure 1.1: The geomagnetic pulsation spectrum 


for the lower frequencies (Campbell, 1966). 
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family of continuous pulsations, have completely 
different origin. They are generated in the inter- 
action of energetic particles with magnetospheric 
Plasma. Pz 2 pulsations: are eleetee to originate 
Peeconance of Alfven waves in the cavity of the 
Magnetospnere in the region of the beginning of the 
meutral sheet (Raspopov, 1968). 

The Pc type pulsations have been found to be more 
prevalent in the daytime hours ( Jacobs 
and Sinno, 1960).  Pc2 and Pc3 pulsations. usually. have 
an amplitude of about 0.5 gamma. Campbell (1967) 
believes that the apparent duration of the signals is 
related to the sensitivity of the recording device. 
They are always present but sometimes too small to 
detect. Pc4 pulsations have an amplitude that ranges 
from 5 to 20 gammas. in the higher Latitudes and ‘an 
event can last from 10 minutes to several hours. 
These are more active in the morning hours with 
maximum occuring around 0800 hours local time. Pe5 
pulsations have the largest amplitude and are some- 
times called og "si Or Giant pulsations, © These 
pulsations occur most frequently near 0600 hours 
and- 1800 hours local time. Pi 2 pulsations have been 
found to, predominate an the mighttime and have a 
Maximum around local midnight. They have also been 
detected in the daytime (Volker, 1968). Pi 1 pulsations 
often appear at £Ehe. beginning of Pi 2. events and also 


appear as riders on ‘the Neng period Pi 2 s.. 
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1.4 Geology and Geophysics of the Area 

ihe ated Imvestigated during this: Study «covers 
Cney central: part of the Province of Alberta, bordered 
by the Rocky Mountains in the west and the Alberta- 
Saskatchewan provincial border in the east, and lies 
between he latetudes 52°N\.and=55°N.e. Geological ly 
this area 1s a part of western Canada sedimentary 
basin. The subsurface geology of the basin has been 
Studied extensively by many petroleum geologists and 
geophysicists with the help of well logs and geophysical 
data. A detailed account of the geology is given by 
McCrossan and Glaister (1966). A brief review of 
geological and geophysical knowledge of the area and 
Posy adjacent Yaneas) wsiciven, here. 

Geologically western Canada may be divided 
PHeLor three distinct wanits., ) (ll) the Canadien “Shacidy 
eGomposed predominantly of ‘crystalline “and meta= 
morpuze mocks: of overall intermediate fcomposi tion 
with local areas of sedimentary and basic volcanic 
rocks. (2)) ‘the Interior Plains uncer lain oy gently 
dipping sedimentary strata of Phanerozoic age and 
is) the Cordillera, a region Of complextly sfolded 
and faulted sedimentary volcanic and metamorphic 
rocks of Precambrian and younger age Ser ei 


intruded by Mesozoic and Cenozoic plutonic rocks. 
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Figure 1.2 shows “the different aimpoortant 
Structural elements and the area covered by the magneto- 
Gelluric sounding. Relative uplifts in the Precambrian 
belt extending across the epicontinental basis in the 


NOTrui-~cast~directron are: also reported by Sikabonyt 


(1957). Faults in the basement rocks are known an 
some areas (Haites, 1960; Burwash, 1965); and in others 
where faults have affected younger rocks. Faults of 


probable basement origin in northern Alberta strike 
northwest vand northeast. Garland and Bower (2959) 
have shown basement faults in northeastern Alberta, 
on. the basis of aeromagnetic data. These faults occur 
with main orientations, northwest and northeast, which 
seem parallel to the known direction of faulting in 
the adjacent shield areas of Alberta and Saskatchewan. 
It has been found that the basement rocks in 
central Alberta are an extension of the Canadian 
Shield rather than of less metamorphosed predominantly 
sedimentary wocks of Precambrian age In tie Rocky 
Mountains. Burwash (1957) has published a 


comprehensive report on the Precambrian based on core 
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Figure 1.3: Lithology of the Precambrian basement 
(Garland¥vand Burwash, 1959), Precambrian 
basement contours and magnetotelluric sounding 
sites in Central Alberta. (The contours 
represent depth to the basement from Sea 


level). 
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13 
samples collected from all Over central Alberta and 
adjoining areas. Garland and Burwash (1959) have 
MmeerDreted- thie Major part of the Bouguer gravity 
enomaly sover, central .~Albertasimn terms; of lithe wdgqiical 
changes in the Precambrian basement beneath the 
sedimentary section. Their lithological map of the 
basement based on gravity data and petrological and 
physical properties of the core samples is shown in 
Fig. 1.3, on which are superimposed the Precambrian 
basement contours and the magnetotelluric sounding 
Sices.. The geographical coordinates jf the magneto- 


pellurrc Sounding.sites are gayengin Hable Zz: 
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TAB IE 2 


ol ae A ee De ee ed ee 
Magnetotelluric 


Sounding Site Latitude Longitude 

ite Alder Flats (AFS) D2 5 TaN dea A) AWN 
24 Bonnyville (BVL) a4% 24)" oN PLO Fe Ww. 
3 Camp Creek (CMC) Bia la aN D4" 477% WwW 
4, Carrot, Creek “(CRC) isi sS orn Lisle Ww 
5. Donalda (DLD) 52° 41' N 112° 37! W 
On Foot ae (PTH) S372 OMY EN LL6>)- 48" W 
es Fox Creek (FXC) 54° 26” N eLCee Se FW 
SaueGreca, bake (GGi,) 53° 30' N 117° 50' w 
9. Meanook (MNK) _ 54° 37! N Pier 20-4 
10. Nordegg (NRG) Serco IN Le 05-9 Ww 
Pipe sNsr Crossing (NSR) Sip to N 6° 40" Ww 
12. Observatory (OBS) Sores. aN abe 42 a Tay 
13. Red Deer (RDR) Ba ON Ia eye a 
14. Smoky Lake (SML) 64°07 aN Le oe 
153 Two aaais (TWH) 53° 44' N is Se." iW 
16. Vermilion (VMN) ee ee oa 110° 44' W 


(Letters in paranthesis represent abbreviations used in 
this thesis). 

The sedimentary rocks, overlying the Precambrian 
basement have a total thickness varying from less than 4,000 


ft. ay the northeass corner Of the warea to greater than 15,000 


ft. in the southwest corner (Fig. 1.3). The thicknesses 
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Figuretl.4: Crustal structuregunder Central 


Alberta (Basham 1967). @ and 6 
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and the electrical properties of the individual layers 
Py ee Scdimetitary “column dirfer from sounding corte 

EO ScOunaing site. “Ihe foeothnills belt in the western 
Sl2ce Gt the “area ts haghly 16lded™and “thrus rau lted, 
tye traults) dipping: generally. to the west and “striking 
NOLUT *Or “northwest. 

Richard and Walker (1959) and Cumming and 
Kanasewrch ~(1966) “have "conducted seismic refraction 
Surveys mn tie p>lains Cl *Soutnern Alberta to study the 
Sart S Crust ifm “this’ region.” ‘Cumming “and Kanasewrch 
estimated the thickness of crust as 45 km. and found 
evidence for a three layer crust beneath the sediments, 
WithsvVelocttres “Gel, “6.5 -and’’/s#2-emesecr* "Basham (1967) 
extrapolated the Cumming and Kanasewich's crustal seismic 
model to the central region (Rocky Mountain House - Alder 
Flats area). The extrapolated model is shown in Fig. 1.4. 
A low velocity channel in the upper mantle in south- 
Central Alberta -is at a depth of 120 Km... and is 
approximately 30 km. in thickness (Kanasewich, 1966). 

1.5 Outline of the thesis 

In Chapter 2, electromagnetic theory applied to an 
isotropic and an anisotropic layered half-space is 
studied. A solution is obtained for the magnetotelluric 
effect in the case Of multi-layered ‘earth, each layer 
containing a dipping anisotropy, and a numerical example 
is presented. Ap short review on the effect of lateral 


inhomogeneities in magnetotelluric sounding is also 
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Given in this. chapter. 

inieChapten >. a-brief description of whe magneto= 
telluric field recording system and the techniques of 
handling and analyses of the magnetotelluric data are 
given. 

Chapter 4 contains magnetotelluric sounding. results 
FOr Sixteen sites in central Alberta with a brief 
description of the results (Or Gach iste. 

in Chapter 5, .an attempt.vs made) to anterpret 
magwerotelluric results fhom central Alberta in “terms 
Of local and regional geology. The effect of the 
source field is discussed. This chapter also contains 


summary ' and suggestions for future work. 
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CHAPTER 2 
MAGNETOTELLURIC THEORY 
2.1 Basic Theory 
The basic theory of the magnetotelluric method involves 
SOLvVing Maxwell's equations with appropriate boundary 
conditions. Maxwell's equations in a source free medium 


can be written as 


oH 


Se San VASE Ze =y la 
ys Bs OE 

Vox, A= J + -€ TE 2 = Lb 

Yom: =>20 Dele =) ee 

V Ei =4. 0 2. legeld 

J =o B, Zoe j= Le 


~ 


where the current density vector J, the electrical 
Conductivity o, the magnetic permeability w, Ehe per-— 


MIECLUVILY «6, Ehe electric: and Magner nertoral weveld 


vectors E and H, are in rationalized MKS units. 


Baking thefcurl of (2.l-la) and substituting: from 
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W2iElS15) “for curl -H: 
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Gin Bar bict aking etheteumh ter (2.1-1b) and substituting 
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for curls 
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and equations (2.l1-lc) and (2.1-1d), the wave equation 


Gan be written as 


I 
Ss 


v2 E 
He ot au, ope 
and 


oH 9-H 


2 S gues = 
Net 3 | le ese = Sao 0. 


Assuming a time variation of the form, 
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Vote Kon = 10, 
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2.1-3a 


Pe AES) ©) 


v-l and w is angular frequency, equations 


Nips: 


The general solution of the wave equation, in cartesian 


co-ordinate system, takes the form 


H H 
= 4 exp (Ker), 
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E E 
O 
~ 
where k.= (k.., Kye Kis 
r= eran sts Pace abe 


and Ho and Eo are constants. 
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In. the frequency range (0.001 Hz - 1 Hz), which is 


of interest in the magnetotelluric method, the dis- 
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placement current term is negligable when compared to 
Eemuctr on current terms (1.e., 6 uw «<< Pp wo) 


Pomcartn materials, and (2,.1—5) reduces to 


eet oO. Ox 2.1-7 
The wave amplitude in a conducting medium decays to 
Wye in a distance -t/Real (k) where this-—distance is 


Known as 'skin depth' and is given by 


6 = / 2 
WWO. 2.1-8 


A skin depth chart for various periods and conductivities 
aseshown in Pigure 2.1. | 
weer vagnevovellunie fields over an Tsotropic;,, Norizone 

boily, layered earth 

To construct a.mathematical model for the magneto- 
telluric effect in the case of a multi-layered 
isotropic earth, consider a cartesian co-ordinate system 
with the following conditions. 

by > The origin: of the -co=ondinare system coincides 
with the surface of the earth, and the z-axis extends 
vertically downwards. 

2). The Garth is divided into (nl) Layers wath 
the cai ten layer extending downward to DHLINVey. 

ay X conductivity value oa and a thickness he is 
assigned to each of the (n+l) layers, where m is the 


layer number. 
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i) tne permeabilvty- and Ehe permittivity « are 
constants within each of the (ntl) layers, and are 
equal to the permeability and permittivity of free 
space, uw and om respectively. 

) LMe Source 1s a current ps Paedbiparnits  -a ata 
height AG from thé=surface of the earth. The geometry 
Cimeae Mode. 1s) depicted, un Big. 2.2. 

The boundary conditions to be satisfied are 1) the 
eengentiral magnetic'-and electric fields are continuous at the 
Horizontal interfaces except: at the current sheet .2) 
at the current sheet the horizontal components of the 
electric tield are continuous but the magnetic: field 
is discontinuous; and 3) the fields must vanish at 
infinity in the z-direction. The boundary conditions 


are expressed as follows, 
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evurrent densaty,in the sheet, Hed Hy and Hy? os are the 
magnetic and electric fields above and below the current 
epect, respectively. 

For simplicity, letathe jeurrentain the™sheet. £4 ow 
in the y-direction with a surface current density Ses 
amps/meter. A current sheet of this type generates a 
beene wave parallel to the y=z plane- with a magnetic 


vector normal to the plane of incidence. With these 


assumptions the wave equation in the mee layer can be 


Written 
ony, s 
Le yt se ~(Cey 2.2 aL 
ae km 


mnie. general solution of which takes the form 
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Since the source is assumed to be uniform 
kK é r ee (Ce) yee, 


ean@.equatvon (2,2-2) weduces to 
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where A_ and Ba are’ comphex constants which can be 

m 
evaluated by the proper application of the boundary 
conditions. The corresponding equation for the electric 


Pieldescan be written, Using W2.1-20) ana (2ea-le), 
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Bpplying the boundary ‘condition, that the tangential 
Magnetic and,éelectric fields are continueus at Z=2 


ere gets: from, (2.2=-3) and. (2.2-4) 
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Pangentral electric and magnetic fields: at each 


Dorlrzontal interface, the constants A. and B. can be 


1 1 
written, 
fay any, 
2. eyes Pee zo - 
A, = pe A, + peee 2, expr i( kjhy) 2.2-7a 
ry ey 
tere ie? 
= a ee I Se See eee et k me 
By 5 A, + 5 °2 | exprit 12) BE 2-2 
At the air earth interface (z=0), where the measurements 


are made, the magnetic and electric field vectors can 


be expressed in terms of A, and Boy 
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o* = apparent resistivity in ohm-meters, 


T = period of magnetotelluric disturbance in seconds, 
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Ba HY = electric and magnetic field disturbances in 
rationalized MKS units 
Z = intrinsic impedance in ohms 


magnetic permeability of the medium. 


= 
I 


The intrinsic impedance at the surface of a 
multi-layered earth was obtained by Wait (1962), who 


defined the impedance in the lowest, eae” layer, as 


nA+1 PPR eis 


and wrote the impedance of any intermediate layer 


Z aE speedtre ( teerrihio (ie vn |) 
geteugtuk m+1 mm r ae censor 
m m m oe + Gta tan - -) ee 
eal 
where Pn = ee 
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Matching the impedance at each boundary, the impedance, 
at the surface, in this case takes the form 


jfiem - Z,+0,k, tanh (kh) 
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2n38i Magnetotelluricifields* over™a’ mul ti= layered! 


anisotropic earth. 

The scalar: impedance as defined in the previous 
section, allows interpretations only for homogeneous 
and isotropic layered earth models. In this section, 
the theory developed in section 2.2 will be extended 


to an anisotropic layered earth model in which one or 
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more layers.,are anisotropic. in, the horizontal (x-y) 
plane. Anisotropy in the horizontal plane may be due to 
banded structure of the rocks of which the medium is 
composed and/or due to two dimensional structures 
which allow the current to flow more freely in one 
direction pers the, other, In such. cases the electric 
field in a particular direction may depend on magnetic 
field variations parallel to, as well as perpendicular 
wou is -durectiion. 

monet (1960) developed a technique to determine 
the effect of anisotropy and/or two dimensional 
structures on magnetotelluric fields using the tensor 


relationships between the magnetic and electric fields, 


H M4 E 

x xX XY x 

H y ¥ BE 

y yx hed yY 2 = 
where Uys = lel vx and vy are complex admittance 


elements. Bostick and Smith (1962) presented a 
computational technique for admittance elements. 
Pokityanski (1961) developed a method to determine 
azimuth “ands coefficient Of anisotropy using) ténsor 
impedances rather than admittances. The tensor 
relationship between the magnetic and electric 


fields in terms of the impedance elements is 
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where Boe Bos are the principal impedances and cer 
ot the additional impedances caused by anisotropy 
and/or inhomogeneity. Madden and Nelson (1964), Swift 
(1967) and Bostick and Sims (1969) have studied the 
properties of these impedance elements and presented 
COmpuLational techniques for their computation 

from the measured electromagnetic fields. Mann (1965) 
has analyzed the problem of a plane wave incident 

on an anisotropic half-space and showed that there 
exist two skin depths corresponding to two principal 
directions of anisotropy. O'Brien and Morrison (1967) 
end=Praus “and Petr (1969) developed theory for 
MULtilavyered anisotropic halt =space: and) presented 
numerical examples for single layer anisotropy in a 
multi-layered earth. Rankin and Reddy (1970) studied 
the polarization properties of magnetoteliuric 

fields in the. case of single layer anisotropy and 
presented a technique to interpret magnetotelluric 
data from the polarization Characteristics of the 
Surface magnetotelluric: fields. “A brier derivation of 


the formulae for impedance elements in the case of 


multi-layered anisotropies is given below. 
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Figure 2.3: Relationship of conductivity anisotropy 


to coordinate axes. 
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Using similar assumptions to those in Section 2.2 
and assuming further that the current flows in an 
folerary direction in the current sheet, and that 
eachelayer is anisotropic in the horizontal plane 
Withee wrinci pal axes x"'.and y', and that the 
major conductivity axis is oriented at an angle 
9 to the x-axis (Fig.)2.3), the wave equations in 
the anisotropic layer, in the two principal directions 


of anisotropy can be written 
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where Ho x and ae m 2re the magnetic field components 
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in the two principal directions of anisotropy in the 


mth layer, and 
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oO and o are the conductivities in the two 
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principal directions x' and y' respectively. The 
general solutions of equations 2.3-1 are 
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2,m 
2.3-3b 
and. the corresponding. electric field components are 


then given by 
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The magnetic and electric field components in the 


measuring. directions, x and y are given in terms of 


E as 
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where oF is the direction of the major axis of 
conductivity (azimuth), as measured in clockwise 
direction from x-axis. These magnetic and electric 
fields can be extended to the surface (z=0) satisfying 
the boundary conditions at each horizontal interface 


and the surface components can be written as 
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The coefficients ay i a Sil and ee Ae defined in 
Appendix I, contain the unknowns Pe peal and a eth 


which depend on the magnetic field polarization in 


the (n+l) layer. These coefficients, A and 
1,n+l 
A. nt] Can be evaluated by the technique described 
v 
below. 


Assuming the magnetic field in the incident plane 
wave to be elliptically polarized, its components, 27H 
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H = G exp (i a) 2515-61 


AO iy Ly 2.3-6b 


where a 1s the phase difference in degrees and G is 
the amplitude recio. From! 2 73-5a, 2.3-5b,.-2.3-0a: and 


was-Ob, using the matrix [R], deétined in Appendix I 


G exp (ia) = (R} cos 0,-R sin ne ar (R, Cos oy = 
R, sin 6,) Ao a+ 2.3-7a 
l= (R) Sin 85 28 R,cos 65) ia (Rosin 0, + R, COs 05) 
are ak 22307 
SS eee So alte eon 
“oh igaeimelar Meare agg CH dah 
where Liye Toye % 43 are the elements of the matrix 
Melee othus A: and A can be obtained. 
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If one attempts to define scalar apparent 
resistivities in the two measuring directions, one 


can write 
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which are not invariant and indeed depend on the 


polarization parameters of the incoming wave. O'Brien 
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and Morrison (1967) computed apparent resistivities 
and developed an interpretation based on the presumed 
knowledge of the polarization characteristics, which 
in ee vary with local time, season and place apart 
Erom local uconductivity structures’ in the earth. 

An interpretation independent of source polarization 
characteristics can be obtained by solving tensor 
fPelation 2.3-2, which can be written out 
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PorsOrven2,s-Li, ler Ol and G be the polarization 


parameters of the incident magnetic field at one 
instant and O.- and G. the parameters at some 
other instant. Then 2.3-1l1 can be written for these 


two instances 
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The four tensor apparent resistivities can then be 


obtained from the four tensor impedances. 
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For a theoretical anisotropic layered earth Zoe sa 
and hence ae = 02 é 
XX VY 
The properties of these tensor elements as 
Summarized by Swift (1967) are as follows. 
FPOrvan rsotropiecyor layered earth 
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when the earth layers contain anisotropy or lateral 


inhomogeneity of two dimensional type 
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For a rotation in cartesian coordinate system the 


rotated impedance tensor must satisfy the relation 
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-sin 6 cos 
gt is *the transpose of fi"and 4 is the angle to whitch 


the axes are rotated in .the Clockwise direction. 


Equation (2.3-15) can be written in expanded form 
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When the measuring axes are aligned with the axes of 


anisotropy 
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and for any other orientation of measuring axes 
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The magnitude of one off-diagonal element is always 


greater than the magnitude of diagonal elements, and 
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invariant under rotational transformation. 

HOE ah anisotropic halt—space or for a ‘single 
layer anisotropy in a multi-layered medium, the 

een | Sees 


direction that maximizes {|Z } or minimizes 


xX yx 


the quantity ae Wa + ge defines the direction of 
the principal axes, or in the case of two dimensional 
Structures the direction of the strike with 90° 
uncertainty. 

In the case of multi-layered anisotropies, where 
the principal axes or conductivity are oriented in 
different directions in each layer, the direction 
obtained by the above criteria is defined as the 
Vapparent azimuth” or “apparent Strike vdivection® 
with 90° uncertainty. This direction is a weighted 
average direction depending on the electrical and 
geometrical parameters of the layers. 

Numerical example for two-layer anisotropy in a 
multi-layered medium 


A simple four-layer model is used to demonstrate 


the effect of anisotropy on the tensor impedances at 
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Figure 2.4a: Tensor impedance elements at the 
surface of a multi-layered anisotropic earth 


in the measuring directionsy(x. andgy). 
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Figure 2.4b: Tensor impedance elements at the 
surface of a multi-layered anisotropic 
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the-esurtace,.. The, top jand bottom, layers .of the earth model 
(Fig. 2-4a) are anisotropic with the principal axes 
oriented in different directions in each layer. In 

Fig. 2.4a the absolute values of the four tensor 
impedances are plotted as a function of period at the 
Surface, in the measuring directions (x and y directions). 
As expected the diagonal impedance elements 2Z and 


xX 


a9 are identical and small compared to off-diagonal 
elements. 

The absolute values of the tensor impedances in 
the apparent principal directions (x' and y') and 
the apparent. azimuth are presented in Fig. 2.4b. 


In the apparent principal directions, the diagonal 


and Z_,., reached their minimum values 


elements Zits! y'y 


(note the scale change for Zt! and Bray le) 

Fig. 2.4b). The maximum separation is observed 
between the two off-diagonal elements in the apparent 
princspali directions... For, short periods the apparent 
azimuth approaches to the azimuth in the second layer 
but for longer periods the apparent azimuth approaches 
30° -athichsi.c..ess than. azimuthiin the thard Laver, 
A5° pe This seffect: is due, to the fact) chat for longer 
periods both individual anisotropies are effective. 
The transition region between 20° and 30° in 


apparent azimuth is marked by the maximum in the 


diagonal impedances. The presence of diagonal elements 


a 


“eben ates aud: eo arovel nod tod bas gos edt a 


wees Lsqioniza vari) ia Ee agouti ling’ ans the 

at. xeyel ee tri - enossosaib dnaxettlb aE) ! 

toe07e4 “dot sift to asiiny 6 sulbeds’ acte = as 

Sat 32s boiweg to aAoitoOnUt B BB: pedo le atts 2eon 

. (erteitvserib y bas x») encisoses PAs itaeberk 303 te é 


5) Pag 1 1 


bos 18 aenetele sons bsqas [pao seep rr pasag 


Leampeto-Ttio at betsqgmos ‘Tame Sits Lepiiariat 


ai #apishbsqmi rosie? sid. 7 eoytev edutoeds | 
bia ("yy Bas *x) enol sser kb Lagtoniag snezee 
yh 


,aa.5 . pit os bedi] 3 ast: ops , ium tip dnote 


Lésops Lb i yENOL 708% t8 tnatoaixg, RSS 


» @enlav mofinta tos tisas berioaax 7 SB 
a | ifn | | 'y? ys a ssa 
th | tok by! Pia bas, eh 102 apnadts etene 
\ | We is 
BR ass Ph "= Bieegiacio ek noi zenages, magn ody + (hs fm 
' ay! 


snexsaes eas rhb sinomele Leregeitb-2id ows ad, he 
| - thexeqas ans eboiseg Hrode 10% -setoesvenip. 
eer bosses. sda nt debut att, a) “asnsaonage 
Bey merisensaae dais, diet & sage ont ‘ebb tisq pay 
| | Wayet by int ortd at drome ‘nade aaah ek 
aegnol ie tess dog and, OF) tiie ek ile 


avis SX6, eeigotioeia evbiviba 


vite nh, 


at oe hrs “OS needed 


end ni momixsm Sits caanate 
met Lanopeth ko voriseonaleat 


oe a. , @ 
7 a Patt | ee ryt 


42 


even ror very=long peraods is due to ‘the fact that 

the apparent azimuth is significantly different from 

the azimuth in the second layer, for these periods. 

2.4 Magnetotelluric effect of dipping anisotropies 
Dipping anisotropies are encountered in sedimentary 

Strata which contain unconfirmities. The anisotropy 

may be due tothe microstructure of the rock. 

Sedimentary rocks for example, are generally more 

resistant in the direction normal to the bedding 

plane than parallel to the bedding plane. When large 

volumes are involved, there may also be present another 

kind of anisotropy called 'macroanisotropy', which 

is caused by a succession of beds with varying 

resistivities which appear to have higher 

resistivity normal to the bedding plane. For 

example, consider a sedimentary column with a number 

of horizontal beds, each with its own characteristic 

resistivity Pin and thickness na For a current 

flowing vertically through a unit cross-sectional 

column, the transverse resistance Ry is the sum 


of resistances met in each of the individual layers. 


R SL NenNH Orh.§ 
m=1 


The average transverse resistivity, which is seen 


by the current flowing vertically through the column 
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Riy m n 
CAS arr A Pn’ y hi 2.4-1 
m=1 m=1 
where H is the total thickness of the layers. For 


current flowing laterally through the column, which 
is the case in magnetotelluric method when applied 
to horizontal layering, the longitudinal conductance, 


which is reciprocal of resistance is 


224-2 
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Unless the resistivity of the individual layers are 
exactly the same, the longitudinal resistivity is 
smaller than the transverse resistivity. This 
dependence of. résistivity on the “direction of current 
flow constitutes anisotropy, and the coefficient of 


anisotropy for layered sequence of rocks is defined as 


—A = YOO, 2.4-3 
This kind of anisotropy may be important in studying 
dipping sedimentary strata by the magnetotelluric 


method. 


This problem can be solved by considering a layered 
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Relationship of axes of anisotropy 


to the measuring axes in the me? layer. 
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half-space in which each layer has a dipping anisotropy 
es *snown in Fig. 2.5. Ohm's law for the anisotropic 


medium shown in Fig. 2.5 takes the form 


Jin 04 0 0 Eon 
Jom = 0 05 0 aaa 2.4-4 
Jon 0 0 oe Eon 


An incident plane wave can be resolved into 
components in mutually orthogonal directions with 
magnetic field vector oriented along (H-polarization) 
and normal (E-polarization) to the strike direction 
(x"-axis) . Since an arbitrary disturbance can be 
resolved into Fourier components of angular frequency 
Owl an exp (1. Wt), time dependence, the electro- 
magnetic field relationship can be written in 
rationalized MKS units (neglecting displacement 


Gurrents),,for the two polarization cases as follows: 


E-polarization 
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where these equations are expressed directly in the 


singly primed system of coordinates. 
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H-pOlarization 


OB WW JE W 
a FO). UE = Mat oh ~ 
x oY 4 2.4-6a 
oH W 
e) E ane 2 
2) yt oz 2.4-6b 
OAL wn 
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where O11 Fy and O, are ENSMEONGUCTIVitTLeS, inexe) “y! 
and z" directions respectively. Expressing these in 


the doubly primed system results in an independent 
set of equations. The solutions can later be trans- 
formed into the singly primed system as will be done 
below. 

3 Since the medium chosen is free from lateral 
inhomogeneities, Ey must be constant with respect 
to y-direction and hence no vertical component, Ho 
of the magnetic field exists. The equations for 
B-polarization are similar to the equations for an 
isotropic, Romedensous earth with conductivity Oj. 


H and Et being thus the solutions of the 
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Helmholtz equation: 
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Hi and Eos im the: am layer are written as the 


solutions of Helmholtz equations in that layer, 
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H = -k - - 
Ve, Mh ole coos il, m2 ERs vg Pan Sab (Se epee 
2.4-9 
3m 
= . Sik a 
ee che setete (ees i,m! 2-1! oa oe an 
U 
WA cena ee 2.4-10 
where 
, ee LOE ICT i Oe 
A Le ws oh jpand 
ie and ues are the complex coefficients. 


Solutions for the electric and magnetic fields in 
the case of H-polarization are complicated by the fact 
that they depend on J, as well as O23. Chetaev (1960) 
solved the H-polarization case for a half-space. Sinha 
(1969) extended Chetaev's work to a two-layer model, 
in which only the bettom layer contained a dipping 
anisotropy. Both Chetaev and Sinha studied the case 
in which no anisotropy exists in the dipping plane. 

In this thesis, a more general model is considered, 
consisting. of a multilayered medium, each layer 
containing a dipping anisotropy, with the strike of 
the anisotropy oriented at an arbitrary angle to the 
oe ee directions, but with one of the principal 


directions, e.g. x', with its associated resistivity 


fe) lying in the horizontal plane. 
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pance there is no gradient ‘of current density-in the 


Zar rectson, by using Ohm's Law (244-4), the above 


equation reduces to 
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From (2.4-6) and (2.4-12) the Helmholtz equation in 


the mo layer for the H-polarization case, can be 


written: 
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The general solution of (2.4-13) in the far layer is 
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From the geometry an Fag. 2..4 
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Equations 
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2 4im9), 


(254-10), 


(2.4-14) 


2.4=20 


(2.4=20) 


specify the electric and magnetic field components in 


the mth 


layer, in the strike and dip directions of 


the anisotropy. 


If the strike of the anisotropy is oriented at an 


angle @, to the measuring axis, x, the magnetic and 


electric field components in the measuring directions 


are written, using (2.3-4), as 
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= m q 
E = i p cos 6. + — = P Sin 6 
ao pli Cina 3), m 2,m 4,m m 
ie 72 2 
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Equations: (2.4-2la) to’ (2.4=21d) are samiitar to 
equations in 2.3-5, hence the magnetic and electric 
field components at the surface of the earth can be 
computed by using the same iterative technique given 
in Appendix I and in Section 2.3, by matching the 


electric and magnetic fields at each horizontal 
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Figure 2.6a: Tensor apparent resistivity curves 
for a three layer anisotropic model in the 
measuring directions for various values 


OL U's 


Suk 


spuodes ‘ GOIYAd 


Ol z0l 


€ Chaat 
l 00:2 '.0-9 ‘0-9 'WYyoodl= I: d 


JIGVINVA =D OE =@ 


E Z : 
BBO eo eee 


4 


Cal 
ya 10:9 ‘0-9 ‘wyol: d= 
(et 
H1d3q 


© 
ssajyou WYO " ALIAILSISAY 


N 


Ol 


LNdaavVddV 


iy. 
ee 


a 
: 
A 


——_ a 1a 

| a 

7 ; | 5 | | : 
A | 8 . 
7 | | ba 
= Ss i < f , 
F *K=n/*O= FS f= = a. < = i ri 
; O=n/70 a) 9 A sm QO! s a aa Se paanareeS* 4 RS 
a mOh= o, eri ay “a Sabo | 7 S 
N VIAAIRAV+D {DE =8 a 3B 
r A eae -_—— : Oe i > 

"0-5 *0=8 Sa = EET SE =X = 
i : 
te. for 


De 


Figure 2.6b: Tensor apparent resistivity curves 
for a three layer anisotropic model (shown 
inePig. 2.6a) dim the strike and.dip directions 


for various values of a. 
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interface. The tensor impedance elements and the 
apparent resistivities can also be obtained by the 
same formulae given in section 2.3. 
Numerical Example 

A simple three layer model is used to demonstrate 
mene effect of dipping anisotropy on the -magnetotelluric 
apparent resistivity. Both the upper and lower 
layers contain a dipping anisotropy. A set of tensor 
apparent resistivity curves in the measuring directions 
corresponding to dip angles varied in steps 0° to 
SO), 2s shown in Fag. 2.6a. ~The model and its 
parameters are also shown in the same diagram. Fig. 
2.66 shows another set of curves for the same dip 
angles computed in the strike and dip directions of 
the intermediate layer. Figure 2.7a illustrates the 
model and a set of tensor apparent resistivity curves 
in measuring directions for a bixedtdip angle a = 40°), 
but varying xX = P3/P>- Figure 2.7b shows the tensor 
apparent resistivity curves for the same model in 
Strike and dip directions of the anieotropy of the 
intermediate layer. In both the models: the strike 
Of the anisotropy is fixed at 30° Benene measuring 


axis, x, and the resistivity along the strike is 
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Figure 2.7a: Tensor apparent resistivity curves 
for a three layer anisotropic model in the 


measuring directions for various values of 
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Figure 2.7b: Tensor apparent resistivity curves 
for a three layer anisotropic model (shown 
i Fig, 2.7a)tim, the Strike and avo uirections 


for various values. of 4. 


B98 


¢. 


2oVins Whitcie esd. Sais: 
. f 


- i! " 
nWone) 2S DOM (5 LLC? Ed. f 


atotissihb git (bas soli 


spuodes ' GOIyId 


suajaw -WYO * ALIAILSISSY LNIYVddV 


56 


When, a.=.0, the results correspond to that of 
horizontal anisotropy with resistivity Py and P5 
along x' and y' directions respectively. The 


resistivity P along.z. -directionshaseno effect 


Since;the,flow of.currentpuis,horazontades »pAnsimilar 
Situation arises when a = 90°, the apparent 
resistivity .1S independent of Pseq When Xkenl 


Gines P5=P3) the apparent resistivity curves are 
independent of dip angle a. The apparent resistivity 
in the strike direction is independent of resis- 
tivities P5 and P3 and dip angle cd. (Figs. 2 56b and 
2.7b). The apparent resistivities in measuring 
directions ,depend.on all.three)resistivatnzesoin 


addi 0n.,.to, strike. and dipsangles sof pbheaanusotsopy . 
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2.5" Magnetotelluric effect in the case of lateral 
inhomogeneities 

Solutions of Maxwell's equations in the case 
of lateral inhomogeneities have been obtained 
for certain two-dimensional geometries. The 
Condition for the, validity of two dimensional 
solutions is that the extension of electrical 
Supuceuresasiould be large compared to a skin 
depth for the frequencies of interest. Further 
the source is considered to be uniform over a 
Similar distance along the strike of the inhomo- 
geneity. 

The electromagnetic field vectors can be 
resolved into components, such that the electric 
field component and the associated currents are 
bare -lel sco .(B-polarization) and perpendicular 
tO, (igo larization), the strike: direction, 

In a cartesian coordinate system, with the 
X-axis parallel to the strike, y-axis perpen- 
dicular to the strike and the z-axis vertically 
downwards, Maxwell's equations for the two cases 
can be written in MKS units, assuming a fea Soke 


time dependency, as follows. 
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E-polarization 
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The scalar Helmholtz equation to be solved then 


reduces to the form 


07°F Lp ul L: 
ees + oes 1WYUOF 0 eS 
where F = E, Lom Hapolarizgation and fF = Hy fOr 


H=podka Giez arbison?. 
Analytical solutions ain the, Hjpolarization case of 
equation }(2:. 5n3)n, fom some siimp.le geomet ri ess.such 
as faults and dykes were obtained by d'Erceville 
and Kunetz (1962), and Rankin (1962). Weaver 
(1963), Blake and Swift (1967), and others 
obtained solutions for E-polarization which were 


incorrect due to improper boundary conditions 
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Se 
used at the air-earth interface. The source 
problem for the H-polarization case is Simple due 
poecne Antanite conductivity contrast’ across the 
air-earth boundary. Thus 
al dee. HOF, ES art, 2° =O 


Brom *(2°.5-'2 ) 


or 
Hoy) = constant. at z/= <0). 

At the surface of the earth the magnetic field 

component parallel to the strike is constant 

with respect to the y-direction even in the presence 

of lateral inhomogeneities. For the case of 

E-polarization the air layer must be included in the Calculation: 

ane tora solution, to be found, this Javer 

must be thick enough’ so that the field in the 

region of the source is relatively uninfluenced 

by the inhomogeneity. 
Neves (1957) first applied finite 

difference technigues to solve some two-dimensional 

problems but unfortunately his solutions suffered 

from errors associated with incorrect boundary 

conditions. In recent years numerical techniques 

have been developed by Madden and Swift (1969), 


Wright (1969), Patrick and Bostick 11969) ana 
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Be nve= and Price (1970a) bd), to Study the effect 
of lateral inhomogeneities of two-dimensional 
type on magnetotelluric fields. Madden and 
Swift (1969) and Wright (1969) used the trans- 
mission line analogy to solve the magnetelluric 
problem, whereas«Patrick and Bostick (1969), 
and Jones and Price (1970a) used direct finite 
difference techniques. A numerical example is 
presented below to show the effect of a two 
dimensional conductivity structure within the 
earth on the apparent resistivities and the 
magnetotelluric fields at the surface. 
Numerical example for lateral inhomogeneity 
This example is taken from Patrick and 
Bostick (1969). The two-dimensional model 
consists of a conductive prism buried in the 
biret layer of a two-layered ecaken. ine 
dimensions of the prism are 1 km x 1 km, the 
thickness otf the. first-Llayer isml. 25m, and 
the top of the prism is 0.25 km from the surface 
Of the earth. The first=lLayer resistryigy wis 
10 ohm-meters, the prism and the semi-infinite 
second layer have a resistivity of 1 ohm-meter. 
The apparent resistivity curves for three 
different sounding sites and the model are 


presented in Fig. 2.8. For the sounding site 
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Figure 2.8: Apparent resistivity curves for H and 
E-polarizations over a two-dimensional 
structure for three different sounding sites 
Situated along a traverse normal to the 


strike (after Patrick and Bostick,“ 19699 . 
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above the middle of the prism, the apparent 
resistivity curves for H-polarization and E- 
POvarization approach the same: Wimi't for Short 
periods. As the period increases ae about 

1 sec. the curves diverge and remain separated 
Over the entire period range above this point. 
Miewsolid curve in Fig. 2.8 is thé apparent 
resistivity curve for the one-dimensional model 
which would be obtained if the prism-were 
infinitely wide. The apparent resistivity curve 
for E-polarization and the one-dimensional curve 
approach the same: limat for, long periods. Fora 
sounding site located at the edge of the prism 
(sounding Site #2), and 0.5 km away trom”the 
edge of the prism (sounding site #3), the 
apparent resistivity curves show a constant 
off-set at longer periods but to a lesser extent. 
As the sounding site moved away from the anomaly 
the apparent resistivity for H-polarization 
becomes larger than for E-polarization. 

The field components computed at Pie er race 
of the earth show that the magnetic field is 
Constant for’ H-polarization, (igs, 2.9), » The 
Glectric field for B-=polarization varies slowly 
across the strike while H-polarization varies 


rapidly... This, implies that for 4°ser.> period, 
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Figure 2.9: Magnetic and electric field variations 
along a traverse normal to the strike for the 
model showh ib Fig./2.8 for 4 second period. 
TE corresponds, to B-polarization case and 
TM corresponds to H-polarizationycase. 


(attcer Patric and Bostick, 196972 
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for which the field components are computed, 
the anomalies in the apparent resistivity are 
reflected in magnetic field for E-polarization 
pucenum the. electric’ tieldt for H-polarization. 
Due to symmetry H, vanishes along the center of 
the anomaly. The peak He occurs close to the 


contacts. 


64 


eae Sts eaonognay bined oi 
orn 9 eae sams tay te 


1h 


abitest log -h wt Sloe? yer iments aids 


| paid rated 4 Abe mies rans? aa 7 
to isunes an? phots soitainke rl veserm 09 | 


eieiee ont sels mcicrob ft tsaq gat .. Vhismnt a 


/ 


65 


CHAPTER 3 
DETECTION AND ANALYSES OF THE MAGNETOTELLURIC 
SIGNALS 
Sek Detection 

Variations in the ‘telluric and geomagnetic field 
are measured at each station using mobile magneto- 
telluric Signal detection system designed and built 
mestie University of Alberta, (VoOzoff et al... 1964) . 

Tae telluric fLeld variations ane derived trom 
potential differences between two copper electrodes 
driven into the ground. The electrodes are spaced 
gees 50 feet) apart an each. Of two, orthogonal 
directions, generally geographic north-south and 
geographic east-west. A 100 wf capacitor is used 
Lo_bdock the D.C. potentials due to, thetpolari zation 
at the electrodes. 

The sensors for magnetic variations are 
induction Coils. Three coils arewburred an mutually 
erenogonal directions to measure Whe Variations 
if) the vertical (Ho), the geographic north-south 
(HY) and the geographic east-west a components 
Of the magnetic field.. Each «oil “consists jor 32,000 
turns of #28GA enamelled copper wire on a 60 inches 


by 0% 75 inch high pérmeability core. 
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recording system. 
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The signals from the-coils and the. electrodes 
are fed into a preamplifier and post-amplifier 
eyecem. Ihe post-amplifiers have discrete gain 
meen e rom 20 1cb co 20nd in steps vot l0rdas Lo 
permit recording Of various signal, levels. The 
amplified signals are then recorded on- magnetic 
tape in analog form using a seven channel Precision 
inscteuments analog FM tape. recorder with a recording 
speed of 15/16 inch/second. WWVB time signals were 
also recorded. A block diagram of the magneto- 
telluric recording system is shown in Figure 3.1. 
The amplitude and phase responses of the magnetic 
@oc telluric systems aré shown in Piguce 3:2. 


ee foto Process ing 


Records approximately three hours in length 
and free from discontinuities were selected for 
digital conversion. The analog signals were 
alzes filtered and sampled by an Analeq to Dagiital 
Converter at a digitizing rate Of a & sanpies (per 
second per chanhel. The analog) tape te layed back 
Se eaesneca Of 30 Mnches/second. Wliierg. ea 
digitizing interval of 32/15 seconds. Digitized 
data sam ll bit brnary form 2m blocks Sr 7200 -ewo 
character words were recorded on two seven channel 


asynchronous, incremental, Kennedy digital tape 
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Figure 3.2: Calibrated amplitude and phase 
response for the telluric and magnetic field 


detection systems. 
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recorders. The data from the two seven track tapes 

were read and stored on nine track tape in IBM 360/67 

Fortran compatable format. 

peeeaePesitral Computations 
Assuming that the magnetotelluric signals are 

samples of stationary random process, the auto- and 

cross-power spectral density estimates for the five 

Components. .wH ».H..» Hoe Eide ail Ey of the magneto- 

telluric fields were computed using the Fast 

Fourier Transform and the technique described in 

Appendix II. The following notation will be used 

to represent auto and cross power spectral density 

estimates and the coherency. 

If X and Y are two time sequences produced by 
measurements on stationary random process, then 
PX(w) = <X*(w) X(w)> = time averaged auto power 

Spectral density estimate. of view series. for 

angular frequency w, 

PY(w) = <¥*(w) Y(w)> = time averaged auto power 
spectral density eStimate.ol Yiseries. for 
angular frequency w, 

PXY(w) = <X*.(w) Y(w)> = time averaged cross power 
spectral density estimate of series X and Y for 


angular frequency w, 


PRY (a) 


GSathy oe (4i)) a FSi ee = coherency between 


® and Y Series for angular frequency w, 
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symbol w of the spectral densities will be omitted 
onthe remainder of this work. 

Digital time sequences containing 4096 points 
Poricns tang roi t3 ,000 sdata tpean ts vand) 10.9 69?zweroBitwere 
weed ifon mos tot ithe spectral «conputatirons-omhe 
maximum lag number was 512. Digital phase shift 
corrections and the removal of instrument phase 
and amplitude responses were carried out in the 
frequency domain. The smoothed spectral density 
estimates were obtained by applying a Parzen 
window in the lag domain. 

3.4 Analyses of Magnetotelluric Signals 

Using the definitions of the preceeding section, 

the Cagniard-type resistivities (scalar resistivities) 


can be expressed in terms of autopower spectral 


density estimates of Ee ne EY and Hoy as 
ere 2 PE 
x 27 U ee 3.4-1 
PE 
Bie a ze 
aaa Stoo 


Phase differences of ay and Booty can be obtained 
from the cross power spectral density estimates, 
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Asfat has been pointed out in Chapter 2, the 
s@alar apparent resistivities are valid only under 
Venyeeim ved, COndLtioOns: for the interpretation of 
magnetotelluric) data. For this reason tensor 
impedances and tensor resistivities were computed. 

owitt. (PfOoo7), Morrrsom etscal., (1968)" and 
Sims and Bostick (1969) have presented methods to 
Solve sche: impedance relation (2.3-—2) for “the 
measured magnetotelluric fields. Consider for 


example the equation 


oe xxix ss oeyily 345 
fromiwhich four distinct equations can besobtained 
involving they various spectraleestimates= in 
terms of the auto and cross power spectral estimates, 


these are 
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of equations, there are six ways to solve for oe: 


Snows 6. olmaher ly. Z eu. 52, can be estimated in 
xy yx MY. 
Six different ways. For example, the six estimates 
One, iz, are 
en 


Pe eas wiliges seeks et, 2. 
haat Aa ys x _Y 
XxX 1 HRT DR Eel Siege rae aest oot cut Tag Ng (ie Sh Sh 3.4-7Ja 
x xX yY y by es 2 yi 
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x x Y x Pe 
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Of these six possible estimates, two tend to be 
relatively unstable, for the one dimensional case, 
Particularly when the incident LLrelds are unpolanized, 
For ‘ehis ‘case ae PELH aut and EERE tend 
toward 2ero, soithe equations (3.4- 76). and Ws wad) 
are indeterminant. The other four estimates are 
Stable and correctly predice Br provided the 


inezdent fields are not too highly polarized. 
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The four stable estimates of each impedance 
element were computed using the are and auto 
power spectral density estimates and an arithmetic 
mean Of the four estimates of each tensor element 
was obtained. These mean values were used in 
subsequent computations. 

The mean impedance elements are rotated 
Shrougn 9.0; tir steps, @t,.5? tostind «che princi ke 
Gvsechions, of ;,conductivity.; » The. direction sthat 
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Strike direction with 90° uncertainty, since the 
geologic strike may coincide either with the 


hasimum, conductivity. or,wicth minimum conductinv ity, 


For the case of a one or two-dimensional 
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must be small compared to unity In iorder that ‘the 
data fit a two-dimensional model. Sims and 

Bostick (1969) used impedance ellipses as a 
eriterion for distinguishing between” ewo- and 
three-dimensional magnetotelluric structures. Each 
impedance element traces an ellipse in the, complex 


plane as the measuring axes are rotated and for 
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two-dimensional models the ellipses degenerate to 
Sireionc lines, which in;the case of the diagonal 
elements Zs and ae Dass through the oragqin:. 

The quality measure generally used for 
magnetotelluric field data over one-dimensional 
Structures is the coherency between orthogonal 
electric and magnetic field components... However, 
in the case of anisotropic and/or two-dimensional 
erouctiures 7, the electriestield component j1 any 
eonesdirection (other than in the strike direction) 
is affected by magnetic field components both 
perpendicular and parallel. Hence the coherencies 
between orthogonal components of the electric and 
Magnetic fields may not be a true indication of 
data reliability. in such cases. “predicted 
coherency" (Swift, 1967). can be wsed as . measure 
Of ther duality of tie deta. , Tie @oredicted 
coherency is the coherency between the measured 
E-field component and the E-field component 
predicted fromy the’ mMagnéeticyfield componentst and 
these averages valhuesof thes impedancetied ementse..c The 
predicted coherencies for both the components of 
the electric field»can besexpressedeas 
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when the predicted coherencies are high, the signal 
to noise ratio is high. 

3.5 Polarization Parameters of the Magnetotelluric 
Fields 

Fowler et al., (1967), Paulson (1968) and 
Rankin and Kurtz (1970) have studied the polarization 
characteristics of micropulsations using auto and 
cross spectral density estimates. The technique 
used by these authors is applied here, to obtain 
the polarization parameters (viz. the degree of 
polarization, the angle of polarization and the 
ellipticity) of the magnetotelluric fields. 

For quasi-monochromatic signals, the degree 
of polarization (ratio of the polarized intensity 


to the total intensity) is given by 


Ro ep 3.4-10 
where J = PX PY - PXY PYX and the X and Y series 

refer:to the components of either the electric or 
magnetic fields. For completely polarized signal 


J>0, R+l and the mutual coherency between X and Y 


series must be 1. For completely unpolarized 
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Signal the coherency between X and Y series must be 
zero. 

tive angle Or polarization Of “the faedd  Cthe 
angie through which. the coordinate axes must be 
Botared' tO align with the major and minor axes 
Oiethne ellipse) weagiven by 


2 Re PXY 
(PX - PY) 


ay Fed 2h 


For positive © values the angle is measured clock- 
wise from x-axis. 

The ellipticity of the field ellipse (the 
bacho Of Minor’ to major axis of the ellipse) is 
given by 

Wem etan 5 sin? is 

((PX-PY) *+ 4 PXY PYX) 
S242 

The elliptacrty ranges “Lromgmytor te circularly. 
polarized signal;* to zero "for’a Pinearly* polanezed 
eignal. Lf © 1s poSitive the Sense ot polarezarron 
is counterclockwise when looking into the propa=- 
Gating wave and the rotation.sénse_is_ clockwise if 
Wy is negative. 


e256 Vertical Magneto Preld, He 


The presence of vertical magnetic fietd 
(assuming plane wave incidence), implies that the 


underlying medium is laterally inhomogeneous, 
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tiene ds; the vertical, magnetic field at a given 
Prequency is) directly proportional to the! gradient 
Po) wene direction normal to the strike of the 
eleceric tield component parallel to therstrike:, 
Hence, a maximum coherency is expected between He 
ancdgcne electric field component parallel to the 
strike, or the direction that gives maximum 
coherency’ between the electric field component 
and Le gives the strike of a two dimensional 
Conductivity structure. 

The coherencies between the vertical magnetic 
tveld and the four horizontal ‘components of the 
magnetotelluric field were computed for various 


rotation angles using the expressions 


PH. E COS?) Go PRE Sin 6 
J, (a ZY 


@lesayls aogier 
ae PH [PE cos*6@ + PE sin*68 + Re(PE E )sin 20) 
Z x y > 


324-138 
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zw ax Zo 


(Oe eel hel ee ; 
oe. PH [PE sin’?9 + PE cos*6@ - Re(PE E_)sin26] 
Z, me 4 x Y 


2. A= 135 
(Similar expressions canbe wrlttem for veoh HH, ,and 
con HOH jj where @ is the» rotation angle measured 
in clockwise direction from the measuring axis x, 
and EL and a are the electric field components 


along the rotated axes. 
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The amplitude ratios between the vertical and 
horizontal components of the magnetic field are 
also obtained for each rotation angle, to determine 
the relative amplitude of the vertical field which 
depends on the conductivity contrast across the 


facerakh discontinuity . 
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CHAPTER 4 
MAGNETOTELLURIC SOUNDING RESULTS FOR 
CENTRAL ALBERTA 

UpEns Chapter |the magnetotel luric results 
(erecive 16 sounding sites, lasted in Pable 2 of 
Chapter il, are presented an alphabetical. order. 
HOreeach station, unless otherwise noted, will be 
Shown the curves and tables as follows: 

a) curves of rotated tensor apparent resis- 
Eivities in principal directions; 
b)/- curves of azimuth, the direction of minimum 
resistivity measured positive clockwise from north; 
c) polarization angles of the electric and 
magnetic fields; 
d) table of values of the coherencies H -Es 
Ge ee Heo Hoy bo and the ratios Ooh and Pe 
in the measuring directions and for the direction 
in which one of the electric field components 
exhibits maximum coherency with He Aeromagnetic 
maps are also presented where available. 

Whe selecti om Criteria: Used §to Jecre tae 
tensor apparent resistivities and the azimuths 
were: 


a)) the predicted coherencies Tor tnexnelec ric 


field components greater than 0.95, 


b) the Skew values tess than or equal to @.2, 
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amnesangles of polarization presented are for the 
degree Of polarization greater than 0.5 and the 
ee eroericity’ less than;0.5. “The second selection 
Criteria for tensor apparent resistivities was 
pelaxed for the Donalda station for which the 
data was slightly-lower in quality and also for 
Nordegg and Observatory where the high skew 
probably indicates a more complex geology. 

The tensor apparent resistivities were 
computed, by the methods outlined in Chapter 3, 
for time sequences for which the original data 
produced high coherencies between EH and 


y 
is ~HY- Coherencies greater than 0.9 for the 


y 

magnetotelluric components generally resulted in 
high predicted coherencies for ES and Be on 
well as low skew values and also smooth tensor 
apparent resistivity curves. For good magnetotelluric 
data the coherencies Apa and Slee were generally 
Tower (<0.8) than the coherencies 

Eh oes ole ed 2 c ONS Da bat High coherency between E, and 

s reaching 1ts theoretical limitvor sunrey 
occurred for the stations in mountains and in 
foothills, which were characterized, in most cases, 


by low coherency for the orthogonal magnetotelluric 


components as well as low coherency between Hy 


d H 
an y 


‘ i I ; De’ ' i : Ae Ao ny ann rei ; ae 7 


any): wot = a8) ‘beanses 


f ert bas SP ciate | 
noLtosl se bre ue eer lea srt a  eenem 
28 esidiyiteies te tbqge! toured 108 sixg 
ie. 7 ae Pn 
: ‘ Ay va 7 
Re eh RN Ue at 
SHt Wot fiw) ros norséie shiemea) ont xo? ‘ge € 
q , i Ne b 
tot oels bas Leon pt vot idepi te as 
wore mpi sty sisaw erodavesed® Snape 
| ha at 
\. Ypolosp xelgncs stom, a aedsotbni 
Siow esidivivelaey oa enbaqs xoeaed on 
wag, } | p+ or 
,€ testqedd ar pent ble o eborsem eds ya , basug 
‘sab Lea eee any Hospi sar ‘epomeupee omis 
i eA te eee | ben ai 
ee naews a0 eetonsxesion ftpid beaub 
a nm | a . " 
| rictt fen O. paris ae 84D. as foassesdo9, _ 
a, 
ff hed fuze 5% vblexonge 4 poe shognos 9B J Leseden 
i Pe AW ip indie 
2.8 6 ya Yo) abasnssasoy besoibes 
: Soa i * 7 a i : 
ae ry. Nd a ee ee a 
toened doom dels bas ean lev, * W 
| ae OF le ae i Dae: 2 
ofiy Cis 3048AvER BoOR a ebeala v 
| | a , Peay. A 
72 4 nae: ¥. } ih, ; wig 
vi Lisweney sxsw ce a brs H>. Hi es fone: 
Pe Oa a Pe, yee ee 
sianersdeo PSs ean enn a 
ian, Peace Dee & 
bas a neewsed Yonoxbries dein 
9 i ar eve ay, 7 
ere | Sent eS 
YER 20 gure! teotds ran 


Jab brs entedavonl 3 enol 
. if Dy ae =r 7 

' _ ; / id i = a, ee i 
: e 428269. Seon fb hes itedoe at 3 
- oe a 

| ee 
, ¢ < a ; a vs | 
at | wll eaenetsd oantes 


eo 


81 

The deqnues..of polarization, angle of polari- 
zation and the ellipticity of magnetic and 
telluric. fields were, computed usdng, the techniques 
described an. Chaptes..3 .«, Jmathe,polerizationsplots 
the number, of, data. sets) msed.for thes analyses; is 
indicated by the number following the station 
name. .An attempt _was).made.to- obtain: a uniform 
distribution of the polarization angles over the 
day. 

Aeromagnetic maps (published by the Geo- 
logical Survey of Canada) are available for some 
of the sounding sites. These depict contours 
Of total magneticwiLeldsyintensttyig Thesélught 
altitude was 1000 ft. above the ground level. 
No correction is made for regional variation, 
which increases ,between 0,..5y.to 3y/mileecirom snorth 
to south and 2.8)y/mile from west to east. 

The anomalies expressed by the magnetic 
contours are generally believed to be dependent 
On. Variations, in. .the magnetic amtensaitiues; of the 
basement rocks, which are, for the most pare due to 
Changes in composition of the basement rocks. While 
strong anomalies are probably due to an increased 
magnetic content in the rocks, weak anomalies may 


well be due to changes in the basement topography. 
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pee Loer tats: (fig 4.1) 


The tensor apparent resistivities show 
anisotropic behaviour for the periods greater 
than 10 seconds. The curves seem to converge 
for periods less than 10 Seconds and approach 
asymptotically to the resistivity of the sedi- 
mentary strata, which is over 4 km thick at this 


LOGcation. The. 0%. values show much scatter when 


a 
yx 


coefficient less than 0.2 indicates a good 


compared to the o0 values. The skewness 
approximation to a one or two dimensional 
Structure. The azimuth of minimum resistivity 
tsivabout +50°. Theyvtelluric f£reld is* polarazed 
in the north-westerly direction. 

Since there 1S ‘no*sign#itcantsvertiecal 
component at this station poor coherencies are 
obtained between He and the other four components 
of the magnetotelluric freid=™ However, Vcoh ey 
and coh HOH are Significantly @greater than 


Coh H E and coh H_H “for N40"? #etation. 
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Figure 4.1: Magnetotelluric results for Alder Flats 
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ZPepeponnyville (Fig. 432) 

Bonnyville is the closest station to the 
Canadian Shield. The thickness of the sedi— 
Memcary Seciion 12s about 125 km. The votated 
tensor apparent resistivity curves show aniso- 
EGopic behaviour for periods Greater than 10 
seconds. The maximum difference between the two 
Curves occurs between 30 seconds and 100 seconds. 
the curves seem to converge for shorter periods. 
Relatively smooth resistivity Ry ele are 
obtained in both the minimum and maximum 
directions of the resistivity. 

The-azimuth for the minimum apparent resis— 
tivity is about +35° for periods greater than 
20 seconds and tends towards zero for shorter 
periods. The telluric, field is polarized about 
N55°W direction... aThe,;magqnetic nana telluric 
fields are non-orthogonal, for most sor the 
spectrum except atqshertygpesiod sends sihere lis 
much. Variation ingthemanaqles of spolarization of 
the magneti¢vfieldeforeperszods lescethen 30 
seconds. 

Mm sLoniii cane Hy component is observed in 
some sections of the data, particularly during 
times of high magnetic activity. «coh HEL 
coh HOES and coh HH, are high in the measuring 


directions. However, the coh HE, OLOps to 
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Figure 4.2: Magneto@elluric results for Wormyvilile 
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minimum for the rotation angle N30°E and coh HUE 

and coh HOH, continue to show a maximum indicating 

a possible strike direction of N60°W. The amplitude 
ratio between H. and Hy varies between 0.1 and 

Ono. The high value for: the ratic ae isa 


consequence of low power in the component. 
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Seeeecamp. Creek. (Pig. 423) 

The tensor apparent resistivity curves show 
Emesocropic behavionr for the periods greater than 
40 seconds. The degree of anisotropy is less than 
edat-Or Alder Flats or Bonnyville. The apparent 
resistivity estimates are relatively smooth for 
skew less than 0.1. 

The azimuth for minimum apparent resistivity 
is about +45° for periods greater than 40 seconds 
ena tends. to 0° “or 90° for Shorter periods. The 
telluric and magnetic fields are orthogonal for 
shorter periods indicating a possible isotropic 
nature of the surface layers. The angle of 
polarization of the telluric field is; about N50°W 
which corresponds closely to the direction of 
maximum apparent resistivity. 

The coherencies between the vertical magnetic 
field component and the other four components of 
ehe Nori zontal magnetite and telfunve a veluseare 
low, however, the coh HEY and coh a drops to 


a minimum for a rotation angle Of N40°E. 


— 


oS 


(£.% .p29): gee ome 


1 Ant Tp tt opeatein mena 


WORD ROP eL VIIVESeISSs bitsy RUT > (oarrnad att 


| . ; ne ; ae 
fey ne ft Cie (Le R511 Sts. aacantaeas Me widelaoy i de 


PISLeSy ANaLEaGs TTA POR dabiniedaiae ent | 


inte yea {) 2 ris tc C37 Fat b.! eet) Le aty BAO St" ep fi ssc 
i 

siteseY sSsptons Bom  “Oeewey = wa Oe eboed pas ‘ 

Yot JBAOBOnoto ers eblest Sigetoss Bae sad Med 4 


set eldieeoa si pattsotbar Shea fedora 
a: 


is oT Stays ly) SUSI aie) aS prise 
i z in f 1 o . 

‘i ‘ hed r Ae 1 y 

OOM Swods 2a bend bl Lend oie ae scien ae emo : 


Oo MOLtI94L Bb vette Os viocols Shaogeexies dod 
a ik Ge ie he iad are ixemy . Ae, 
UWIPAM itso LIAS OD) ose eS Lt et eterloo pe! 


fee 


Om BINS TOON OD, iyOT cSt cr sgn ier fF KOEMOD ‘tei 
Ste ablert (ocr l Pad pis Ea he 3 sits || 
4 panne 2 ve 


3 eh ON Gite af ee) Sirs a cu nos ay 4 Tavewod ,wol 


Om 2 ee ; et bis Wy ; a“. eo 
eA" UGH 26 alone OOF EC ye BOD iumen in s 


89 


Figure 4.3: Magnetotelluric results for Camp Creek 
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Tensor apparent resistivities in the 
principal direceone'. 

Azimuth (the direction of minimum tensor 
apparent resistivity). 

Polarization anglesa teri tellur.e.and 
Magnetzc fields. 

Vertical magnetic field analyses’ 


results. 
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4. Carrot Creek (Pig, 434) 

Me Tensen apparent resistivity curves for this 
station are similar)ito the curves for Camp Creek. 
MNReEManisotrop1c effect is pronounced for periods 
greater than 80 seconds, although the presence of 
Busisenoti1ceable fOr periods: greater than 20 
seconds. There is a considerable amount of scatter 
in Pry even for skew values less than 0.1. 

Rotation of tensor impedances gave +50° 
azimuth for minimum apparent resistivity for longer 
periods and tends to 90° for shorter periods. 

There is a considerable amount of scatter in 
azimuth. The angle of polarization of the telluric 
field distributed around N45°W for long periods 

and close to east-west for short periods. Magnetic 
ana telluric faelds*\ ane. orthogonalr forishorter 
periods. There is a considerable variation in the 
angle of polariLzationlom Chel magneties Timeicey 
although it is mainly oriented along the N20°E 
direc ts.0n:. 

Significant levelulot He activity is observed 
er this site. Coherency analyses show £hatecoh Siti 
and coh HOW, are greater than coh HOE. and coh cease 
further <coh HE, and coh HOEY drops to a minimum 
for a rotation angle of N30°E, giving a probable 
strike direction of N60°W. The amplitude ratio 


between H and Hy Varnes. Dotween U..? anc, 0.5 Lor 


tris. rotation angle. 


90 


1 4 VE Lk ey 
ay eer’ i i es 
7 1D :* i. f is 
{ , hv we 
ne ih 
i Db i‘ 
{ i 7 4 
i 
ae 
\) 
A hi 
( * a ‘ 
i 7 i 
of bd he ' ~ * , 5 mar 
Tn tap Lark A it x a I tb 4 “ 


” t 
, » 1 . [ oe 1 - ! 
he 
t i } 1) i er, 
yes em thee . , fy ah ew i 5 Tart i 
SHOU Leu) Lo SOOO VURL Ota ah eee ot ot ee wv 


in ’ 
oi Lie . 
i) {{ <i i Ld Ss SCY SS tae 
Tag al 
i Pei kd 1. oa TO), Be 
; ) : 
; 8 - 
| [AV Wade "fig 


fol tot yeas eles tS 19h age’ en aa tm 


i 


it 


ehorned Yaewode oR OR ie abips Sas- &. 


1 r } Cane ; i. 
Ysiitsoe to snvoms 2 aid: stabi Jaoo se ai ox 

| Fran Vein 9 is ah 
otiviies eAdiie i ae! edt As an: 


ie 


iboitase prol Sot bv Pe Note bedvatatee bf 
\ ny % ‘Ww BY, : : Welw ue as 
sitanpeM  .ebobisa ae ey Fea sabe fond seols 


1 e4 torte 45% ienopod3 40. S38, aie kw te 


Sid fh TO LIE Cavs 
7 ne Mad 


BLS Le bors 


H ; , ¥ or | tu ql * yi 7 a 
“ Ria ; ; ves a ' of cares A 
OSV aeco: 2 £ \wed Le Lip paw BY : i 2 


TH, doo ser, wole 2Ssey Dene rye 


HUH toe) bas | ti-d> ede ae 
ba - bm ; je y i 7 ny) 7 


PILI eee Ot sy ey ae 


? i, 7 


otdadots & preva’ aed 


Figure 4.4; 


oa 


Magnetotelluric results for Carrot Creek 

4 Tensor apparent resistivities in the 
Princ#pal sawrections. 

li. ~Azimuth 4themdirection Of minimum 
tensor apparent resistivity). 

Lii., Polawazattomngamgles for telluric and 
magnetic fields». 

iv. Vertical magneti@ Field analyses’ 


results. 
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The Guality of the data, for Donalda was not 
very good and did not meet some of the editing 
requirements. The skew values for all most all 
of the data sets analysed were greater canon 
although the predicted coherencies were greater 
than 0.95 for some data sets for certain period 
bands. Rotated tensor apparent resistivities for 
skew values less than or equal to 0.3 are smooth 
along the maximum apparent resistivity direction. 
There is a considerable amount of scatter in 
Oxy Values and they continue to drop even for very 
long periods. 

Rotation of tensor impedances gave an azimuth 
between +30° and +50° for periods greater than 
20 seconds and evidence of decreasing for shorter 
mDervods. The angle of polarization of the tellumic 
f,eld as about N40°W and, there 1s much scatter jo 
shorter periods. The degree of polarization of the 
magnetic field was less than 0.5 in most cases. 

Coherencies between the vertical magnetic 
field and the electric field components are low, 
however, good coherency is obtained between H, and 


H for longer periods for some data sets for the 


rotation angle about N40°E. 
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Magnetotelluric results for Donalda 


als Le 


seas 


tas 


Tensor apparent resistivities in 
the* principal direetions. 

Azimuth (the direction of minimum 
tensor apparent resistivity). 
Polarization angles fior teliluric 
and magnetic fields. 

Vertical magnetic field analyses 


results. 
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FOOot Hills ‘station is ‘situated in the Rocky 
Mountain foothills belt. The basement depth is 
unknown in this area. The tensor apparent 
resistivities showed too much scatter and did not 
Satisfy the selection criteria and hence they are 
not presented here. The power level in the magnetic 
and telluric signals was reasonably good for 
several records but the coherency between magnetic 
and telluric fields was low (<0.8). The skew values 
were greater than 0.5 and the predicted coherencies 
were wwess tian 0.9. Rotation of tensor ampedances 
GavemcamevazLmuch pout to0' | but Signiticant 
amplitudes for diagonal impedance elements were 
obtained even for this angle. 

The degree of polarization of the magnetic 
faeidiieor one of the two records chosen for 
polarization analyses is less than 0.5 for most 
of the periods. The electric field 1s linearly 
polarized in the north-westerly direction. 

The coherency analyses with the vertical 
Magnetic field show that there 2s a considevabic 
band of coherencies greater than 0.8, between ae 
and aly for the unrotated data. A rotationsof 
N50°E improved coh H HY at the expense of coh ao 
An anomolously high value for coh HOEY appeared over 
most of the range Of periods an the rotated 


direction. 
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Pigure 4.6: Magnetotelluric jresults for Foot, Hills 


i) Polarization angles” for tie telluric 
and magnetic fields. 
if. Vertical magmetic field anaibyoes 
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7. Fox Creek CPign 4.7) 

The rotated tensor apparent resistivities show 
an aMisotropic behaviour for periods’ greater than 
20 seconds, with Oxy Values scattered for: peraods 
Greater "chan 200 seconds.) The Oxy’ Values Sate scoO 
decrease for periods longer than 200 seconds. For 
Short periods the apparent resistivity curves tend 
to approach 10 ohm-meters in value. The azimuth 
which tends to zero degrees for periods less than 
10 seconds approaches +60° for longer periods. The 
angle. Or polarization "of the tellurie field is 
about N35°W for periods greater than 100 seconds 
but for short periods it 1S oriented along an east- 
West direction. The magnetic fielad@as polarized 
along N40°H for longer periods and about NoO-E ror 
shorter periods. 

Good coherencies between the vertical component 
Of the magnetic. field and the four hormpzoncal 
components of the magnetotelluric field are obtained 
fe the Measuring directions, for Long pemloads.¥ “Vor a 
rotation angle of N50°E, coh HOE, dropped to its 
minimum and coh He and HH reached their maximum 
Values implying a strike direction of N40 WwW. 
coh HOH, is vhighitorvall the periods Tor igs 
rotation angle whereas coh se Starts te decrease 
for periods less than 25 seconds. The amplitude 
ratio between He and a for this rotation angle, 


Varies with period from 0.3 to 0.6. 
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Magnetotelluric nesultss for Fox Creek 


ie Tensom apparengaresistivetwes in 
the principal directions, 

iis Azimuth (the direction of minimum 
tensor apparent resistivity). 

iii. Polaraezgation anglesstcor telluric 
and magnetic fields. 

ys Vertical magnetic field analyses 


results. 
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8. Gregg Lake (ee 457.8) 

The Gregg Lake station is situated within 40 
miles of the central range of the Rocky Mountains which 
Have a north=-westerly strike in this area. The 
Srevaviomvor the station 1S. about 4,000; feee and. 1 a6 
Surrounded by mQquntains up to 6,000 feet in height. 

The coherencies between the orthogonal magnetic 
and telluric field components were low and the scalar. 
resistivities were highly time dependent with much 
scatter. The tensor apparent resistivities had much 
Scatter and did not satisfy the selection criteria 
and hence are not presented here. 

The electric field as linearly polanxezed in the 
north-westerly direction. The angle of polarization 
Gr the magnetic f1leld varied with period, We rs 
approximately N30°E tor long periods and abeue NGS" 
for shorter periods. 

Sida tioantly Greater oF H was observed for 
Snorter pervods than tor longer periods.  Coberency 
analyses show that iy is coherent with Es ve Hy 
and f for shorter periods in the measuring, @irections, 
eEnd’for & rotation angle of N50°E the coherency between 


H -E. and H -H_ dropped to their minimum whereas the 
x Zz 


Gone hae ana t£neccoh aoe showed an increase. 
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Figure 4.8: Magnetotelluric results for Gregg Lake 


i. Polarization angles for telluric and 
magnetic fields. 
Tk Vertical magnetic field analyses 


results. 
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9. Meanook (Paige Ag) 

The Meanook magnetotelluric sounding site is 
Situated within a few hundred feet of the 
Dominion Geomagnetic Observatory. Niblett and 
Sayn-Wittgenstein (1960) and Srivastava and 
Jacobs (1964) conducted magnetotelluric measurements 
and interpreted the scalar resistivities in terms 
of upper mantle conductivity anomalies. The 
tensor apparent resistivity curves from the present 
analyses show an anisotropic behaviour for periods 
Greater than 30 Seconds. The apparent resistivity 
estimates are smooth for skew less than 0.1 for 
long periods, but there is some scatter at the 
short-period end of the spectrum. 

The azimuth of minimum resistivity is found 
to be about +35° "for periods greater than, 30 
seconds and about +15° for shorter periods. There 
ms much scatter in the angle of polarization of 
the magnetic and telluric fields for shorter periods, 
butrtor Long peniods the telluric field 1s > polarized 
about N55°W whereas the magnetic field preferred 
N35°E,. The direction of maximum apparent Tesistivity 
andthe angle of polarization of the telluric field 
are the same for these periods. 

The vertical component of the magnetic field 
is coherent with the four horizontal components of 


magnetotelluric field in the measuring directions, 
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Figure 4.9: Magnetotelluric results for Meanook 


i Tensor apparent resistivities in- the 
principal directions. 

iM. Azimuth (the direction of minimum. 
tensor apparent resistivity). 

iii. Polak zation angles for telluric 
and magnetic fields. 

Lis Vertical magnetic field analyses 
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for long periods. For a rotation angle of N40°E, 
coh HOE, and” coh BH dropped to their minimum 
values and coh ek and coh HH, reached their 
maximum values, which is consistent with a possible 
Strike direction of N55°W. 

An aeromagnetic map which is available for this 
area shows a strong north-east trends a little to 


the east of the Meanook site. 
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10. Nordegg Giants gee i) 

The Nordegg sounding site is in the foothills 
ef the Rocky Motintains at an elevation of 4,500 
eee.) )1t) 1s surrounded by mountains ap to 6,000 feet 
magn on all sides. Thirty miles to the west of the 
sounding site lies wae First lange. of the Rocky. 
Mountains with an elevation up, to 10,000 feet. 

This and other ranges in the Rocky Mountains have a 
MOrcA-westerly trend in this area. Magnetoteliuric 
measurements were made at this site during the 1968 
and 1969 field seasons. The coherencies between the 
orthogonal components of the magnetic and telluric 
f2elds were low (<0.8).for both years" records 7% fhe 
coherency between the orthogonal components of the 
EFelluric tield was very high. “Predicted coherences 
were good only for afooupilesofsmecondss fom both 
these years and the skew values were greater than 
#70. Although the tensoeri appanenitnt es: situ ithy 
Curves show anisotropucibeha viouMigs Pt ts not 
advisable to draw any conclusions from these > 
results, because of haghliskew values. The 

azimuth for minimum resistivity was found to be 

+90° (not presented here). 

The telluric field is linearly polarized about 
N30°wW, for both years, with slight deviation for 
1969 records from 1968 records, which may be 
There 2s mach scatter in the angle 


imercrumental. 


of polarization of the magnetic field but it seems 
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Figure 4.10: Magnetotelluric results for Nordegg 


ule Tensor apparent resistiwmities in 
measuring directions. 

LiL.) ' Polarization anG@les fer Gelluric wnd 
magnetic fields. 

iii. Vertical magnetic field analyses 


results. 
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to be near an east-west direction for most of the 
period bands. 
Coherencies betWeen the vertical component 


of the magnetic field and Eos Bey ies and ny are 


y x 
high for periods greater than 100 seconds and for 
a@ rotation angle of N40°E the coh HOH reached 


their maximum without much change in other 


coherencies. 
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11. NSR Crossing (Fig, 4.11) 


The station North Saskatchewan River Crossing 
(NSR Crossing) is located at the crossing of the 
North Saskatchewan River and the Banff-Jasper 
highway, in the central ranges of the Rocky 
Mountains. The station is surrounded by north- 
westerly striking mountain ridges with peaks up 
Co 1-000. feet. | 

Very low amplitude magnetic signal was recorded 
at this site. The coherencies between the orthogonal 
components of telluric and magnetic fields were 
low (<0.7). The predicted coherencies were also 
low, and hence no reliable apparent resistivity 
curves were obtained for this site. 

The polarization analyses of the magnetic and 
welluric field Show that the telluric freldmic 
linearly polarized in north-westerly direction tor 
all periods whereas the angle of Seen nh OF of 
the magnetic field is about N20°W, for perzeds 
Jonger than 100 seconds, and very scattered for 
shorter periods. 

Very low H, activity was recorded at this 


site. 
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Figure 4.11: Magnetotelluric results for NSR crossing 
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2. Observatory Cigna LZ) 

Magnetotelluric results for the University of 
Alberta Geophysical Observatory were obtained both 
from 1968 and 1969 data. The skew values were 
high for this location. The rotated tensor apparent 
resistivities presented here, are for skew values 
Pescethan Or.equalrro 0.3. The apparent 
resistivity curves show an anisotropic BenGnons for 
periods greater than 8 seconds but seems to converge 
£LOm Shorter periods. The Ov Value contain less 
Scatter than Oxy Values. 

The azimuth for minimum resistivity changes 
Prom +50° for long periods to *90° for shorter 
Periods. The angle of polarization of the telfuric 
field is N40°W for long periods and there is much 
Scatter, for shone periods, but it iS oh between 
north and north-west. The magnetic field at 
NeO°E darectionwtor long periods. (For “snort periods 
it scattered around N35°E. 

The direction of maximum resistivity and the 
angle of polarization of the telluric fiedq 
coincide for longer periods. 

For data sets with reasonably good He Lor 
long periods, the coherency analyses show that H. 
is coherent with the four horizontal components of 
the magnetotelluric field in measuring directions, 
and for a rotation angle of N50°E, coh HEY and 


coh H H. dropped to their minima whereas other 
Za Sy. . 
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Figure 4.12: Magnetotelluric results for Observatory 
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the principal directions. 

Azimuth (the direction of minimum 
tensor apparent resistivity). 
Polarization angles for telluric and 
magnetic fields. 

Vertical magnetic field analyses 


results. 
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coherencies were virtually unaffected. 

The tensor apparent resistivity curves 
presented here are similar to the curves obtained 
in the measuring directions for 1966 data (Rankin 
and Reddy, 1969). The difference in the direction 
of principal axes obtained in both the cases may 
be due to change in the direction of magnetic 
Pieicd pOlarizatrvon, Jas well as less accurate 


analysis techniques applied to 1966 data. 
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H3. Red Deer ( EGGS. Anes 

Magnetotelluric data for the Red Deer sounding 
Site was collected by Dr. W. J. Peeples “as a part or 
his magnetotelluric investigations in Southern 
Alberta. The tensor apparent resistivity curves 
Show an anisotropic behaviour for periods greater 
than 40 seconds. The curves are isotropic for 
periods less than 40 seconds and seem to approach 
10 ohm-meters for shorter periods. The minimum 
apparent resistivity estimates have much scatter 
for periods greater than 200 seconds. 

The azimuth) £6r minimum resistivity is around 
20° for periods eeneee than 40 seconds and varies 
between 0 and 90° for shorter periods. The angle 
Sf polarigation of the electric: field 1s, about 
N50°W for the entire band withtsome ‘scatter-sdt “ene 
short period end. The average angle of polarization 
for magnetic fiéld is ‘around N40°E Pp bute there is 
much scatter. The magnetic and satéeliurive fields yare 
approximately orthogonal over the entire spectrum. 
The direction of maximum apparent resistivity and 
the average direction of ttelluriemtteld tare 
coincident. 

The amplitude of the vertical component was 
low for this station. Coherency analyses show that 
i is more coherent with ay and Hy than with ES. and 


H_, for long periods, and for a rotation of N30°E; 
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BLGULe4.]-3-——Magqnetotellurie scestltsS for Red Deer 


i. Tensor apparent resistivities in 
the principatgdirections. 

a. Azimuth (the direction of minimum 
tensor apparent resistivity). 

iii. Polarization ang¥#es for telluric 
and magnetic fields. 

iV. Vertical magnetic field analyses 


results. 
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14. Smoky Lake (Pa Ggs-4.. 14) 

The tensor apparent resistivity curves in the 
principal directions show anisotropic behaviour for 
(ae entire spectrum with a greater degree of 
BEnaSOcropy at the long period end» The apparent 
resistivity estimates are relatively smooth, 
although Pay have more scatter than Poy Values. 

The curves may converge for very short periods. 

The azimuth of minimum apparent resistivity 
varies from +90° for 10 second period to +30° 
for 300 second period. The angle of polarization 
Of the telluric field is around N50°W for long 
periods and tend towards N25°W for short periods. 
There is much scatter for short periods. The 
average direction for the magnetic field is 
N20°E, but there is a considerable deviation 
BLOM thls direction. 

Coherency analyses of the vertical magnetic 
field with EL ve Hy and show that H, is 
highly coherent with Er Ey and ee for periods 
greater than 50 seconds, in the measuring directions. 
Por a rotation angle of N50°E, coh HEY dropped 
to its minimum whereas coh oe and coh H, 
remained more or less unchanged, implying a possible 


strike direction of N40°W for periods greater than 


50 seconds. 
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Magnetotelluric results for Smoky Lake 
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Vertical magnetic field analyses 


results. 
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Ley 
An aeromagnetic map available for this area 
shows two very strong positive anomalies within 
five miles of the sounding site, one in south-east 
direction and the other in south-west direction. 
A weak east-west trending anomaly is present north 


-of the sounding site. 
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Figure 4.14a: Aeromagnetic map for Smoky Lake area 
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The tensor apparent resistivity curves in the 
principal directions show an anisotropic behaviour 
for periods greater than 30 seconds. The maximum 
apparent resistivity estimates are smooth compared 
with minimum apparent resistivity estimates which 
eisplay considerable amount of scatter for 
periods greater than 100 seconds with lowest 
value around 200 second period band. For shorter 
periods the curves seem to approach asymptotically 
to 5 Ohm-meters, indicating highly conductive 
Surface layers. 

The azimuth for the minimum apparent 
resistivity is about +40° for periods greater 
than 50 seconds and tend toward zero degrees for 
Shore periods. “Thetélectric fleld is polarwzed 
about N50°W direction for long periods, (> 260 
seconds) and is about N40°W for short periods, 

The angle of polarization. of the magnetic fierd 
is around N20°E for long periods and about N40°E 
for short periods. The electric and magnetic 
fields are close to orthogonality for short 
periods. 

A low amplitude vertical magnetic field 1s 
Bbserved for this Station. The coherencies 
between Hk: es, and tea & are greater than 


0.8 for periods greater than 150 seconds, in 
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Pigure 4-153" Magnetotelluric’ results” fometwo-ttits 


ay, Tensor apparent resistivities in 
the prin¢gipal directtoge-. 

Sigs Agimuth (the direction of minimum 
tensor apparent resistivity). 

chp saa Polarization angles for tel¥urie 
and magnetic fields. 

inva Verticalfmagnetic Eield analyses 


results. 
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measuring directions and coh HE, dropped to its 
fimenamum for arrotation angle..of N40°E. 

The aeromagnetic map for this area shows a 
strong north-south lineament five miles west of 


the sounding site. 
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Figure 4.15a: Aeromagnetic map for Two Hills area 
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or. a Vermilion CBG Arde cbr) 

ine, Vernmilvon Sounding sate ts located close 
to “the Alberta-Saskatchewan provincial border. The 
Pabekness Of the sediments at this site ls about 
Peis “The tensor apparent resistivity curves in 
ee principal directions show an anisotropic behaviour 
Hor (erlods greater than’ 20 seconds, and seem to 
Ponverge cto 6 Oohm-—meters -for shorter periods... The 
maximum apparent resistivities are relatively 
smooth compared to the minimum apparent resistivities, 
which have considerable scatter for periods greater 
Gham LOO seconds. 

The azimuth for minimum apparent resistivity 
is about +40° for periods greater than 30 seconds 
ane wend tO Shark. towards +15° ror @shoerter Pemiods.. 
tie angle of pollarizgatroy Of~ane Celie tia 
Varies from N40°W for’ 40 Seconds” period to N6Ss-wW 
for perilods abowt. s00™seconds. "Por “sirovuer perlods 
tne average direction or the teliurtC Pieter us 
about N70°W. The average direction of the magnetic 
field about N20°E, witth* some scatrer stor short 
periods. The average’ direcevions#or magnetic and 
telluric Eields are orthogonal for short perrods.. 

Low amplitude H, is observed at this location. 
Conerency analyses for one record with significant 
Hoy for periods wWwreates than 100) seconds; “snow that 
He is coherent with Ey Ey! Hy and ie in the measuring 
directions and for a rotation angle of .N50°%E the 
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Figure 4.16: Magnéetotelluric resudes for sewni lion 


auch Tensor apparent sgesistivieies in 
the principal directions. 

ant Azimuth (the direetiom of minimum 
tensor apparent resistivity). 

sia pee Polarization angles aor feelluric 
and magnetic fields. 

Wee Vertical magnetic field analyses 


results. 
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The most important characteristics of the results 


presented” in Chapter” 4-are: 

1) The tensor apparent resistivity curves in 
ives principal directions for all the stations in 
the plains of central Alberta exhibit anisotropic 
behaviour. 


2) The degree of anisotropy varies from station 


tor station, with the maximum “occurring at Bonnyville 


(2), a station close to the Canadian shield. The 
degree of anisotropy 1S°a minimum for Carrot 
Creek (4),°Camp Creek: (3), Fox Creek (7) and’ for 
Meanook (9). 

3) Higher than normal skew values in the plains 
are obtained for Donalda (5) and the Observatory 


(ES thee i ie In the foothills and the Rocky Mountains 


anomalously large skew (> 1.0) values are obtained. 


4,eThe tensor apparent resistivity cunves in 
tnetarrection of the principal axis of maximom 
apparent resistivity are smooth comparedgto the 
curves along the principal axis of minimum 
resistivity. 

5) The azimuth (the principal axis of minimum 
apparent resistivity) lies between +30° and +60° 


for periods greater than 30 seconds, ft6r all the 
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stations in the .plains, -The azimuth is inder-— 
Minette, for shorter periods, except for Smoky 
Lake (14), where the anisotropy exists even at 
shorter periods. 

Oy the, electric field.is linearly polarized 
about.a north-westerly direction for periods 
greater than 30 seconds both at mountain stations 
POCO alm Stations... Bon short veriods the 
Peveamizelion.of the electric fiéld is Jess con= 
Solent ~ although 2 tends to lie in theygnorth— 
west.quadrant for most.of.the times for which 
the analyses were made. | 

7) The magnetic field tends to align nbelmete! 
north-north east direction for longer periods. 
For shorter periods there is a considerable 
Verieataon in «the angleuot polarization of sine: 
ellipse. ,The magnetic field was elliptically 
Celarized withoan ellipticity vess, than 05 iA 
MOS Cases’. 

SB) A signiticant ,vertical componen tof sence 
Magnetic,field was ,observed for Jonger periods 
when the magnetic activity was high. Hy tends 
to be more coherent ~withethe electric frela 
component along the major axis of apparent res- 
istivity and with the magnetic field component along 
the minor axis of the apparent resistivity than 


with the other components. This characteristic of 
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He is reversed or absent for stations in the 
foothills and mountains and also for Donalda (5). 
ey, FOr the Stations in the Olains. the tensor 
epparent resistivities in the measuring directions 
(north-south and east-west) exhibited the 
mablest degree Ofpanisotropy, except for Smoky 
bare (14)., where Significant anisotropy was 
Sosenved in measuring directions at short periods. 
Tensor apparent resistivity. curves along the 
major and minor axes of apparent resistivity and 
the directions of maximum apparent resistivity 
for all the twelve stations in the plains are 
presented in Pigs. Sak, 5.2°and 5.3 respectively, 
Fiei-aApparent. Te6StSiivuty cul Veo Inekig., eu 
show a wide spread from the high values at 
Bonnyville (2) to the low values at Two Hills (3). 
This effect may be due to variations in thickness 
Bia/or the nesistivity,;contrasts.) A. set of 
theoretical apparent resistivity curves for a 
two layer earth model with varying thickness 
auc resistivity of the first Layervare presented 
im Fig. 5.4, as a comparison to the observed 


curves presented in Fig. 5.1. 
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Figure 5.1: Tensor apparent resistivity curves along the 


majOk axis, @& anasotnopy for central Albemmta:. 
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Figure 5.2: Tensor apparent resistivity) curvesMalong the 


minor axis Of anisotropy foriicentral Alberta. 
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Figure 5.4: 4 Theoretical apparent resiStivity curves 
for a two layer earth model for various thicknesses 


and resistivities of the first layer? 
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Tensor apparent resistivity curves in the minor 
direction also show much variation (Fa@y S'i2 Weaeon 
onevstation to “another, °'For!some™stations (eg. 
Donalda, Red Deer) the curves continue to drop 
evemvat very Veng* pera eds, “itor other stations 
they increase with period tending to approach 
asymptotically the curves along the major axis. 
Phere-as*much scatter ane the minimum apparent 
resistivity values. For some stations, 
viz., Bonnyville and Two Hills, the apparent 
resistivities reached their lowest at inter- 
mediate periods. 

The scatter in the apparent resistivities 
along the minimum. resistivity axis 2S expected. 
Since the principal axes are chosen such that 
the diagonal terms of the tensor impedance matrix, 
Zeealeand'2o ,- +? ares /a® minimun {zero in thew idea 1 
case). Thus one can write as an approximation, 


Pnethe*prinerpal “ditection® (x"*and y*) 


The maximum resistivity results from a larger 
value.of E. for any. given, value, of, the-orthogonal 


component of H. Since the Glectmic tiel clas 
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polarized with rélakively small eerie remicagteya et aie 
most cases, approximately along direction of maximum 

apperent resistivity and H is bolarized 
approximately orthogonal to E (close to the 
minimum apparent resistivity direceiom)y 
both Ey and Dan are Of low amplitude and hence 
the sagneal to nétse ratio is poor. VThis, results 
bona welarcively large scatter in the apparent 
Cesroeivity Curves along sk” . 

ine velectric and) magnetic tield polarization 
engles for all the sixteen stations are plotted 
ited. Gao And 546 vVespectivel yy. —Didgterent 
symbols are used.to indicate the case of 
eillipreleuty woess, than 0 Ss dndegreatem than Orso. 
There are only a few cases, generally for short 
periods, where the electric field has ellipticity 
Gueater /Eman 0.5. .tThe magnetic field hasean 
average ellipticity larger than Ceat oF pene 
electric field but also has only a@ few ‘cases wien 
ellipticity larger than 0.5. nese plete wineicar. 
that for periods greater than 80 isecondes tie 
tellairicec fveld vs polarized approximate? a wierta. 
West, most of the values lying weengny 2205 oF 
this direction. For shorter periods there 1s a 
great deal more scatter witn the direction 
lying mainly in the northwest quadrant. The 


magnetic field polaniveblon- angie 1s much more 
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Scattered with an average direction approximately 
north-north-east. 

Devs’ generally’ believed that the angle of 
potarization“of the telluric® field is mainly 
Comer ouled by tive hocalVand/or TeqionaL geo= 
e1ecuruce Structures. ine magnetic field 
polarization observed at the surface may be 
controlled by both the generating mechanism 
Of geomagnetic pulsations and the geoelectric 
Swe UCeUre.- ~“A2ybpin ~(1967)", van van and 
Zybin (1967) made an extensive study on the 
effect of laterally inhomogeneous geoelectric 
structures on the magnetic field polarization 
and concluded that the preferred direction of 
polarization in the natural magnetic fields is 
determined by the geoelectric structure of the 
underlying medium and daily and seasonal variations 
iM-polartzation are~connectéd with micropulsation 
mechanisms. The magnetic field component 
paralter*to the strike of the structure is 
unchanged (Section 2.5) whereas the magnetic 
field component normal to the strike is increased, 
resulting a change in the polarization character- 
istics of the total magnetic field, which is 
measured at the surface. A two dimensional 
structure tends to orient the magnetic field 
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ihe Main Cectonmi@ Mndts Of the’ Geqion are 
oriented in a north-west direction parallel to 
the Rocky Mountains. A general cross section 
passing through Edmonton in a north easterly 
GPreceion, is’ presented anv rig. 5.7. The Vocations 
Gr itive o1 the magnetotel luricrsounding sites 
are also marked in the same figure. 

Sedimentary rocks of low resistivity 
overlying a gently dipping highly resistive 
basement, blocked on the west by more resistive 
Precambrian metasediments and Phanerozoic 
sedimentary rocks belonging to the Rocky Mountains, 
eoOnstitute aamajor,.<north-wests+styiking, .two- 
dimensional geoelectruesstrnucture...l[hbe 
resistivities «Gisethe individualslayers aor the 
sedimentary Golumn vary trom a few Ohm=meters wo 
approximately one hundred ohm-meters. The 
resistivities of the basement rocks may very from 
several hundred ohm-meters to several thousand 
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Figuremo.n<s 59 Geological “eross ‘Section-along 
north-east and south-west passing 
through Edmonton, and approximate 
location of magnetotelluric sounding 


Sites . 


git ‘a hy ie . 
2260 teow ~igwoe bab sese~seon ) 


; ' / i ee Ney a ‘ 
otensneiage 4 wes inp mb ‘dpuose 
‘ ee 4 vs sia | 
POLOmi Oe ait’ 1 flatevenpan to eistannesi f 
o iat Pei H . ay jo " iat ee ; i" 4 


\ Neeaie me 


WY) 
a 
<< 
— 
z 
a 
Oo 
= 
> 
6 
2 ©- ae oS Fy 
Ss 1 mere yt! ; BF if 
S =X SST aU ULA TAL Ui { 7 
SS ee: Tse Horperave 
© a o = 
° N wn N wa ~N 
1 


ASVS 
‘VIT1V 


NOLNOWG3 


499} puosnoyy ‘'NOILVAII3 


WwW 

a 
wn U 
SG) 
Om 
FV 
Oz 
.. & 
© 
I= 
ae 
an 
z= 
oO 

fe 
a < 


PRECAMBRIAN SEDIMENTS 
dl AND SEDIMENTS OF LATER AGE 


e . 
° . . 
° ° . 
Or ge 


a17..5 


even 


iio arom con ones ing 
ter ot pabdel 
=e Fs de® 


ed 


clare 
ogee ee feat 


3JADe 


oman 


2UO3M0} ‘ValaaM ADIN. 
24D08 DIHISOMATIM OMA 


 eTMamidze narnmaneae, 
JOA AITAI IO VTMAMIGRe ie 


139 
detailed knowledge is available on the resis- 
tivities of the rocks in the mountains and 
foothills, however, magnetotelluric measurements 
at Nordegg (Fig. 4.10) gave an average resiStLivi ty 
of a few thousand ohm-meters. 

For a two-dimensional structure of the type 
described above, currents tend to flow more 
freely along the strike in the sedimentary 
basin producing a maximum electric field along 
strike for given resistivity - “This may explain 
the north-westerly polarization of electric 
field (Fig. 5.4). The Secondary magnetic field 
produced by induction trends normal to the 
Strike and this secondary field together with 
the primary field produces a preferred direction 
epouL Weo fF for the total Magnetic freide (rigion5). 

The coherency analyses involving the vertical 
component, H, of the magnetic show that He 
is coherent with the electric field along the 
strike direction and with the magnetic field 
normal to the strike, which is consistent with 
the theory Of the magnetotelluric fields for two- 
dimensional structures. The amplitude ratio 
between AS and the magnetic field component 
normal to the strike was less than unity for 
The ambiguous results obtained with 
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to a more complex geoelectric structure which 
also produces high skew values for these Sitacwone. 

The anisotropy observed in the tensor apparent 
resistivity curves at each site may be due to 
1) Che anisetropy anesthe rocks below the sounding 
Site and/or 2) lateral inhomogeneity due to 
structure, such as the one shown in Fig. 5.7. 
Anisotropies of various degrees have been observed 
in sedimentary rocks (Keller, 1968), but the 
@xistence. of such anisotropies in igneous rocks 
iS doubtful. According to Rankin and Reddy (1969) 
a*veststivity contrast of lito 00 talong ‘the 
principal directions would be necessary in the 
basement and sub-basement rocks to explain observed 
anisotropies in central Alberta. It is most 
unlikely that the anisotropy shown by the 
apparent Lresistivity curves 1s due to the tocal 
anisotropy in true resistivity. 

For the second alternative it will be useful 
tourecall! that for most. of the stations aa the 
plains, the skew values (<0.2) are consistent 
with a two-dimensional structure, which will be 
that one that is presented in Fig. 5.7. 


assumed 


In Fig. 5.8 are plotted the directions of 


maximum apparent resistivity for three periods, 


200, 50 and 10 seconds. For the mountain 
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Figure 5.8: Directions of maximum tensor apparent 
Fesistuvaty for the central Alberta 
stations. (dashed arrows represent the 


polaniZatrwonmangles fergretlurive treldyy. 
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Eelluric field iss plotied,;)sinoeune frejaiab ie 
apparent resistivity curves exist. These directions 
are approximately parallel to the strike for 200 
and 50 second periods, and hence the curves of 
maximum and minimum resistivity can be associated 
with the electric.and magnetic field polarizations 
respectively. The two southern most stations at 
Red Deer ;(.13), and Donalda (5). rare probably 
influenced by the Sweetgrass Arch (Fig. 5.6). 
For, che, .k0,..second. period it .can be sseem tthat athe 
results are scattered but the prevailing directions 
are close to the measuring directions which is 
constitent with an isotropic upper section 
(Sedimentary) of the model. 

The two-dimensional example of Patrick and 
Bostick. (1969.). (Section). 5).andsiahseg ther, wom 
of Swi ft..(1967).,~cand Wraght (1970) ,lGhow iia 
near,surface, gecelectric structunes ‘affecti che 
H-polarization apparent resistivity curves even 
at, very long periods »(tens) of; thousands) of 
seconds) .».However ,cing the resultsspresented 
here the effect is limited to a few thousand 
seconds after which the apparent resistivity 
curves for H- and E-polarization tend to converge 


and approach asymptotically to the resistivities 


of boettom.layerss Magnetotelluric measurements 


at longer periods are necessary to validate the 


ahove etatement. Unfortunately no computational 
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results are available for the type of model 
presented in Fig. 5.7, except for some preliminary 
work by Peeples (personal communication) which 
demonstrates that this type of model yields 
apparent resistivities similar to the one 
peesentedian this thes is., 

bateral inhomogeneity in the case of F- 
polarization has less effect on the apparent 
Resistivity curves than in the case for H-= 
polarization. The curves for E-polarization may 
give an approximately valid vertical conductivity 
structure at a given site. Hence, an attempt 
is made to interpret E-polarization curves in 
terms of depth to the basement and resistivities 
of the basement and sedimentary rocks. 

A one-dimensional model was derived for each 
Giese the tour sounding sites, Viz, Alder Flacc ye 
Bonnyville (2), Smoky Lake (14) and Vermilion (16),, 
The depth to the basement is fairly well known for 
all the four sites. Well logs are available 
within a few miles of each of the sounding sites 
and these were used to check the thickness and 
the resistivities of the sedimentary layers. 

There are very thin highly resistive layers in 
the sedimentary section, but a study of various 
one-dimensional models with layers of various 


thickness and resistivities showed that the effect 
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Figure 5.9: One-dimensional models for four 
magnetotelluric sounding sites in 


Centra leyloes tas 
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of thin layers such as are Present in central 
Biberta 1S insignificant in comparison with the 
scatecer ian the data. 

Several one-dimensional models were 
constructed for each site in order to match by 
trial and error to the observed tensor apparent 
resistivities for the E-polarization. Application 
of computer graphics, which has been tried at 
Some; stages Of interpretation, or adopting inverse 
techniques will reduce the labour in obtaining a 
proper model. The one-dimensional models which 
gave curves in good agreement with 
the observed data for all the four stations are 
presented in Fig. 5.9. The model curves and 
the observed data are presented in Fig. 5.10. 

The models yielded a low resistive sedimentary 
Section for Smoky Lake (14), Vermilion (16) and 
Bonmyvyilie (2) and a more highly wesistive 
Sedamentary section for Alder Flats .1) 3 “Whus 
result is consistent with the well log data. The 
low Jace eee for sedimentary columy £0or Smoky. 
Lake, Vermilion and Bonnyville may be due to 
relatively thin or absent Cambrian, Devonian and 
Mississipian formations, in these regions. The 
basement depths are read from well logs and used 


as invariant quantities in constructing models for 


each site The different values for basement resistivities 
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Figure 5.10: One-dimensional model apparent 
resistivity curves and the observed 
apparent resistivity curves for 
E-polarization, for four magnetotelluric 


sounding sites "An -central Alberta. 
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op caurterent: locations may be die to It theological 
variations in the basement. Garland and Burwash 
(1959) predicted from gravity studies that these 
lithological variations extend up to 8 km in the 
basement. 

The magnetotelluric data and the one- 
dimensional models presented here, validates the 
application of the magnetotelluric method for 
studying the basement profile as well as for 
studying local and regional geoelectric structures. 
The magnetotelluric theory for a layered earth 
model show that the results presented in Fig. 
5.9 are not unique and the work on the effect of 
lateral inhomogeneities show that the apparent 
resistivity estimates for E-polarization cases 
will be greater than for one-dimensional models 
(eee the examples in Section 2.5), If the 
magnetotelluric fields are affected to a 
significant extent by the lateral inhomogeneity 
of the type shown in Fig. 5.7, the basement and 
crustal resistivities must be lower than those 


obtained from one-dimensional theory. 
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5.2 cummary 

The magnetotelluric results in central Alberta 
are affected by the geometry of the sedimentary basin 
and adjacent Rocky Mountains. The apparent resistivities 
fOMeshortéer ‘periods, corresponding toa skin depth 
in the sediments, display no obvious anisotropy, 
indicating that the sedimentary section is free from 
both microscopic anisotropies and lateral 
inhomogeneities. 

The varying degree of anisotropy from station 
to station may be due to local variations in basement 
composition and structure. The inconsistent results 
obtained in the Rocky Mountains are probably due to 
complexity in the local geology. An interpretation 
of these results will require further knowledge of 
the effects Of three dimensional structure on 
magnetotelluric fields. 

While a non-plane wave type source may contribute 
to, the vertical component of the magnetic field, Ho 
the high coherencies between H. and the horizontal 
components of the magnetotelluric fields, particularly 
with telluric field components, indicates that H. iS 


a secondary, induced field and the plane wave 


approximation is Wa elCl. 
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For One-dimensional models deduced for the 
tour -stations, shown anFig. 5.9 on the basis 
of E-polarization results, are consistent with 
the geological data, although these models are 
not unigue. No evidence for a highly conducting 
layer correspogding to low velocity layer in the 
upper mantle is obtained from results of this 


Work. 
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2-3 Suggestions for Future Work 

The LrésultsPpotsthiseworktinditaterithati the 
magnetotelluric method has considerable potential 
for studying regional geoelectric structure as 
well as for mapping the basement in a sedimentary 
basin. An extension in the frequency would be 
most valuable for achieving this potential. 

Present day modelling techniques for two- 
dimensional structures are computationally too 
expensive, for trial andi error interpretation. of 
magnetotelluric data. Hence, less expensive 
two-dimensional modelling techniques as well as 
knowledge of the effect of three-dimensional 
geoelectric structures on magnetotelluric 
fields, might make the magnetotelluric method 
an important tool in geophysical exploration. 

The results presented here indicate that 
a-rerudyrot* source fields and the 
effect of complex geoelectrical structures on 
plane and non-plane wave sources are essential 


to the development of interpretation. 
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interface, the coefficient vector [cl], in the 
first: Layerrcan be written as: 
fee Peto eee (Al) 


where [C41! LS; thewcoeki-cient vector limithe 


ene layer and 


[R] = [X 


ee 


Since the (n+l layer is assumed to be of 


intunitce thickness, the coefficients 6 


ease 
and eel 0p) veannG: 
~ 
Res 
0 
ike = 
n+1] = Breas 
ie | 
Expanding (Al), one gets: 
A = xr A eG A 
Leb ee ie a olae 13 2,n+l 
B = enn A ap a A 
aes Da Herp gael 23 PA rap Nae 
A anes A ae A. 25) 
a ae erik 33 nh 


vt) 


“ i a iM? 


Ae 
ey ve 

ey menace 
iv my 


tp 
2 


ae | 
me) ea 


pe mot nw ea a> = seed 


ae) r 


2) ign oe 


siobtaeen: edd seeendoidd 9 


3609. eds ek, Nag h 


oye 8 (ce) ‘ena’ 


4 sy 


9 ow ‘ hte | 


ere ft Pala | has 


Bas, wowed | ts) 


5 


Heat 


@ 4 = 


AD 


A4 


where 


Ay 1 ’ Bay ; Ay y and Bo are the coefficients 
Pieter tre t layer wand TE 1» Yay and 
C43 are the elements of the matrix [R], which 


are dependent on geometrical and electrical 


properties of the.layers. 


i 


30 staem 


i J 
» > iy feats 
5 at ww 


ene bolts 


A5 


APPENDIX TI 
SPECTRAL ANALYSIS 

The formulae used to compute the auto and cross 
power spectral density estimates are summarized 
here. Detailed derivation of these formulae can 
be obtained from Blackman and Tukey (1958), 
Bendat and Piersol (1966), Gentleman and Sande 
(1966), Bingham et al. (1967), and Jenkins and 
Weates: (19-67). 

Consider a sample record x(t) cf finite 
length T, with N, even numbers of equally sampled 
points and A the time interval between two 
SuccesSive points. The transformed aes values 
are then given by 
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by Fourier series representation 
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The coefficients AL and BL are given by 
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The estimated (raw) autocorrelation function at the 


displacement T = rA is defined by the formula 
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Where t = rA, r is the log number and m is the 
maximum log numbers. Ry is the raw estimate of true 


autocorrelation Ry. The autocorrelation function 


may take on negative as well as positive values of 


Tiyocancd 


i = R(t). 
A normalized. value of autocorrelation function is 
Obtained by dividing Rta) by R, (0), where 
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The autocorrelation functions can be computed more 
efficiently from the Fourier transforms using the 


relationship (Gentleman and Sande, 1966) 
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R, (x) = = C(r) 
N ‘i N 

Citi a en ee ay Cm) eae x* (k) 
n=] Nica 
x aes a 


Where xX 1s the Pourher transform of xe) and * implies 
Complex conjugate. 
Power spectral density estimates 

The power spectral density estimates can be 


obtained from the correlation function 
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Since only the positive frequencies are of physical 
interest, the physically realizable spectrum is 


defined as two times S,.(£) for positive frequencies 
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k is the frequency index also called the harmonic 


number. We can also Write “Px as a function of 


frequency index 
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Smoothed power spectral estimates can be obtained 

by multiplying correlation function with appropriate 
weighting function. The smoothed spectral estimates 
are written as 
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is used in this thesis. 
Cross power spectral density estimates 

Consider a second record y(t) of the same 
length as x(t) with same number of sampled points 
and same sampling interval. The sampled cross 
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amplitude spectrum estimate is 
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